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Abstract: Two effects, jet broadening and gluon bremsstrahlung induced by the propagation of a highly 
energetic quark in dense QCD matter, are reconsidered from effective theory point of view. We modify the 
standard Soft Collinear Effective Theory (SCET) Lagrangian to include Glauber modes, which are needed 
to implement the interactions between the medium and the collinear fields. We derive the Feynman 
rules for this Lagrangian and show that it is invariant under soft and collinear gauge transformations. 
We find that the newly constructed theory SCETq recovers exactly the general result for the transverse 
momentum broadening of jets. In the limit where the radiated gluons are significantly less energetic 
than the parent quark, we obtain a jet energy- loss kernel identical to the one discussed in the reaction 
operator approach to parton propagation in matter. In the framework of SCETq we present results for 
the fully-differential bremsstrahlung spectrum for both the incoherent and the Landau-Pomeranchunk- 
Migdal suppressed regimes beyond the soft-gluon approximation. Gauge invariance of the physics results 
is demonstrated explicitly by performing the calculations in both the light-cone and covariant gauges. 
We also show how the process-dependent medium-induced radiative corrections factorize from the jet 
production cross section on the example of the quark jets considered here. 
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1. Introduction 

The start up of the Large Hadron Cohider (LHC) has reinvigorated strong theoretical interest in the physics 
of energetic jets [1, 2]. At present, many of the new developments in the Quantum Chromodynamics (QCD) 
theory of hard jet and particle production are motivated by the need for precise evaluation of the Standard 
Model background to signature processes for new physics [3, 4, 5]. Such accuracy relies on factorization 
theorems [6, 7] that allow to separate the perturbative hard scattering part of the cross section from the 
non-perturbative parton distribution and fragmentation functions. 

Recently, Soft Collinear Effective Theory (SCET) [8, 9, 10, 11] has emerged as a powerful new tool to 
address large processes in lepton-lepton {t^ +^~), lepton-hadron {t^ +p) and hadron-hadron + p or 
p + p) collisions. SCET in conjunction with QCD factorization, which has been proven in this framework 
for a number of processes [12, 13, 14, 15, 16, 17, 18, 19], is particularly suited to improving the precision of 
multi-scale calculations^ through the resummation of large Sudakov type logarithms [17, 20, 21, 22, 23, 24]. 

An important multi-scale problem is presented by the production of jets in reactions with large nuclei, 
such as lepton- nucleus + A), proton- nucleus {p + A) and nucleus- nucleus (A + A) reactions. In these 
processes the energetic quarks and gluons must traverse a region of dense nuclear matter of (9(5 fm) 
and their interactions in the medium induce broadening and a new type of radiative corrections that can 
significantly alter the corresponding jet cross sections and shapes [25, 26, 27, 28]. Preliminary results from 
the Relativistic Heavy Ion Collider (RHIC) and new results from the LHC indicate that these effects may 
indeed be observable in A+A collisions [29, 30, 31, 32, 33, 34]. It will be natural to use SCET to describe 
the collisional and radiative interactions of the jet in the medium when the typical transverse momentum 
of the partons in the jet is comparable to the size of the momentum exchange with the medium and much 
smaller than the jet energy. So far, only broadening of the final-state parton in semi-inclusive deeply 
inelastic scattering (SDIS) [35] has been considered in an effective theory of QCD. Reference [36] argues 
to present a less model-dependent result but its gauge invariance remains to be demonstrated explicitly. 

To set up a general SCET framework that can describe the collisional and radiative processes induced 
by the propagation of an energetic parton in strongly interacting matter, to derive the operators that de- 
scribe the momentum space evolution of the propagating quark or quark+gluon system and to demonstrate 
the gauge invariance of these results is the main goal of this paper. By comparing our findings to previous 
calculations of parton broadening [37, 38] and energy loss [39, 40, 41, 42] in nuclear matter we identify 
areas where the effective theory calculation will be able to improve the accuracy of existing computations. 

Our manuscript is organized as follows: in section 2 we review very briefly the basic concepts of 
Soft Collinear Effective Theory. The kinematics of jet-medium interactions beyond the static scattering 
approximation is discussed in section 3. We identify the regime relevant to high energy jet production 
in hadronic reactions with large nuclei and elucidate the possibility for constructing an effective theory 
to describe parton propagation in dense QCD matter. The gauge-invariant Lagrangian for this effective 
theory is constructed in section 4. Feynman rules are derived in the light-cone and covariant gauges. In 



""^We refer here to energy, momentum or mass scales. 



section 5 we evaluate the transverse momentum broadening of jets, induced by their cohisional interactions 
in the strongly interacting medium. Radiative processes are discussed in section 6. Our focus here is on 
the soft gluon limit, when the energy of the emitted gluon (cj) is much less than the energy (E) of the 
quark that splits: u <^ E, for comparison to previous results. We demonstrate the gauge invariance of the 
jet broadening and energy loss results in section 7. We deduce the kernels that describe the broadening 
and medium-induced bremsstrahlung as a function of the quark interactions in the medium in section 8. 
An application of the reaction operators for collisional and radiative processes is also discussed in this 
section. The extension of radiative energy loss calculation beyond soft gluon approximation is presented in 
section 9. We also show how the process-dependent medium-induced radiative corrections factorize from 
the hard jet production cross section. Our conclusions are given in section 10. We have moved some of the 
background technical discussion to appendices. 



2. A brief overview of SCET 

SCET [8, 9, 10, 11] is an effective theory of QCD which describes the dynamics of highly energetic quarks 
and gluons. The relevant physical scales in this effective theory are the hard scale Eh ^ Et ^ £^cm, the 
jet scale Ej ^ p± that describes the width of the jet in momentum space and the scale of soft radiation 
Eg ^ Aqcd- The power counting parameter of SCET A defines the hierarchy between the hard, jet and 
soft scales. We use the version of SCET, which is sometimes referred to as SCETj, in which the scales are 
Eh - Ej - A^ and Es - A^. The de grees of freedom in SCET are collinear quarks {^n,p)^ collinear gluons 
(An^p) and soft gluons (Ag). Collinear particles have momentum in light-cone coordinates pc ^ [1, A^, A] 
and soft particles ps ^ [A^, A^, A^], where we define our light-cone notation in appendix A. All other fields, 
such as hard quarks and gluons, are integrated out from the QCD Lagrangian. Their effect on the dynamics 
is contained into the Wilson coefficients of the SCET operators, which can be calculated using a standard 
matching of full theory onto effective theory. In order to avoid confusion we note that what we call soft 
gluon mode in this paper ps ^ A^ is sometimes called ultrasoft, while the soft momentum is defined as 
Ps ^ A, see for example [11]. However, below in section 4 when we define the momentum scaling of the 
source, one of our choices corresponds to ps ^ X and we call it a soft source. 

The Lagrangian of SCET [9] arises from substituting into the QCD Lagrangian ij; = e~^^^ t(jn,p and 

integrating out the small component of ^pni where in — ^i^n^ — ^i^n and ipn — ^ ^n- The result 
for the collinear-soft Lagrangian is: 



>CsCET(Cn, An,As)=^,n 



in-D 



^in^CYM{An,As) , (2.1) 



CYM{An^ As) = ^tr { [iD^^ + gA^^^^, iD^^ + gA^^^^] + ^g.f. , (2.2) 
CG.YiRd-f^{[^Ds,.A^^^^]]\ (2.3) 
£G.F.(LCG(6)) = ^tr{6^A^^j' . (2.4) 

Here, the covariant derivative D contains both collinear and soft fields: iD — id ^ g {A^ + Ag), while Ds 
includes only the soft gluons: iDs — id ^ gAs. Thus, the collinear and soft modes are coupled in the 
SCET Lagrangian. In the first term of Eq. (2.2) the summation over label momenta is understood 
implicitly, and in Eq. (2.3), Eq. (2.4) summation over the label momentum q is understood implicitly. We 



have written out explicitly the gauge fixing terms for the covariant and the light-cone gauges. The ghost 
terms are omitted for brevity. 

A key ingredient of the SCET formulation is the BPS transformation [11]. This transformation consti- 
tutes a collinear field redefinition which involves soft Wilson lines and removes the interactions between soft 
and collinear fields in the Lagrangian of SCET up to the power corrections. Such decoupling is essential 
in the proof of factorization theorems in SCET. The BPS transformation redefines the collinear quark and 
gluon fields: 

Cn,p — ^Ci^p 5 (2-5) 
^a,M_3;a6^WM,6^ (2.6) 

where the Y{x) and y^^{x) are Wilson lines built out of the soft fields in the fundamental and adjoint 
representations correspondingly: 

Y{x) = Fexp(^ig J ds n- + x) T^^ , (2.7) 

y''\x) = P exp (^ig J dsn-Al{ns-\-x) (-i/^^^)^ . (2.8) 

To derive the Lagrangian of SCET in terms of the decoupled colhnear fields ^n\An'^ one needs the 
following key identities: 

yr^yt ^ ybarpb ^ ^2.9) 

Y^n-DY ^n-Dn. (2.10) 

In particular, the last equation removes the interactions between the soft gluons and the collinear quarks 
which is contained in the covariant derivative in the left hand side of the Eq. (2.10). As a result, after the 
BPS transformation we obtain: 

^SCET i^n, An, A^) = (^n, ^n) + l^s (^) + l^cs (Cn, ^n, A,) ^ , 4°^) + Ls{A,) . (2.11) 

The couplings between the collinear and soft modes is removed from the SCET Lagrangian. However, 
in order to preserve the gauge invariance one has to put soft and collinear Wilson lines into the external 
SCET operators. The collinear Wilson line is defined in the following way: 



Wn[x)^V^xp[-iQ \ d5n-A^(n5 + x)r^ . (2.12) 



3. Kinematics of the in-medium jet interactions 

In this section we describe the kinematics of jet interactions in QCD matter. Our goal is to identify the 
typical momentum exchanges between the energetic incident partons and the medium, which in turn will 
help us construct an effective theory for these interactions. 

The typical jet transverse energies that are phenomenologically interesting at RHIC are in the range 
Et ^ 10 GeV - 50 GeV. At the LHC this range is extended to Et ^ several hundred GeV. On the other 
hand, the mass of the particles in the medium is on the order of or larger than 1 GeV. In cold nuclear 



matter this is the mass of the nucleon tun = 0.94 GeV and in the quark- gluon plasma this is the mass 
of dressed partons (x fi — gT^l + Nf/6. Binding effects inside nuclear matter can significantly increase 
the effective mass of the recoiling particles. One of the main goals of this section is to study how the 
elastic scattering cross section depends on this mass and to identify the kinematic configurations that give 
a dominant contribution to this cross section. 

The cross section for the elastic scattering of two particles of masses mi, 1712 is given by the standard 
expression: 

- , . \, , ^^'^I' ^o^ia'gp JllL i^^)'sHpi +P2-PS-P4) . (3.1) 

We work in the rest frame of the medium, which for the cases of a quark-gluon plasma at mid-rapidity 
coincides with the laboratory frame. The same is true for fixed target p+A experiments. In this frame the 
flux factor conveniently reduces to Flux = 4m2Pi. Let us denote by 6 the scattering angle of the incident 
parton {6 — Z{pi^ps)) and by ^3 the physical solution arising from the energy constraint mi)+m2 = 

£"3(^85 ^1) + ^4(^4 —Pi—Ps^ ^2)- In this paper we are interested in the case: mi < 1712. In this case there 
is a single physical solution ^3. Substituting in Eq. (3.1), we obtain: 

^ ^ ^ ^1 (3 2) 

dn 647r2 pim2 [psiEi + 7712) - P1E3 cos(6>)] ' / ' ^ ' ^ 

We now discuss the specific channels that may contribute to (|Mp). To study jet broadening and 
energy loss, we will be interested in forward scattering where, on average, the direction of jet propagation 
is not significantly altered per interaction. Such scattering is dominated by the t-channel gluon exchange. 
One can, of course, write down ^/-channel and 5-channel diagrams but these describe hard backward parton 
scattering and isotropisation processes rather than transverse momentum broadening and energy loss. For 
simplicity, we give the specific example of quark-quark scattering: 

2. _ 4 /2\ 2(i/2 + s'^) + 4mf (2t - 2m^ - mf) + 4m^(2t - 2mf - ml) + Smjm^ 



Here, one can conveniently take the mi 0, m2 limits. We note that if the interaction is of finite 
range, i.e. we have an exchange gluon of mass m ^ /i, the only change is t ^ t — fj? in. the denominator of 
Eq. (3.3). 

We can now study the m2 dependence of the average squared matrix element (|Mp), the flux and 
phase space factors (PS, Flux), and the differential cross section da/dQ. We chose a massless incident 
quark of energy 10 GeV. The mass of the recoiling particle varies from 1 GeV to 10^ GeV and we use 
^ = 2 for the coupling constant. The top left and right panels in figure 1 show the differential phase space 
(normalized by the initial flux factor) and matrix element versus cos(0) with the leading m^ dependencies 
taken out. While for cos(0) = 1 the curves come together, the difference for small values of m2 away 
from the forward region can be many orders of magnitude. In contrast, in the differential cross section 
da/dQ this variation largely cancels everywhere except for the backward scattering region. We show this 
result in the bottom panel of figure 1. The insert illustrates the remaining subtle mass dependence in the 
forward scattering region. It is quite remarkable that this residual variation is less than ±50%. At the 
level of the integrated scattering cross section, the differences are even much smaller. For incident partons 
of £; = 10 GeV at RHIC and E = 100 GeV at the LHC we obtain: 
a(m2 = 1000 GeV) - cr(m2 = 1 GeV) , ^ a{m2 = 1000 GeV) - a(m2 = \GeV) . ^ 
a(m2 = 1000 GeV) Irhic-13%, a(m2 = 1000 GeV) Ilhc - 2% , 




cos(e) 



Figure 1: Kinematic factors and mean squared matrix element for the t-channel scattering of interest are shown 
versus cos{0) in the top left and top right panels, respectively. These exhibit strong dependence on the mass m2 of 
the recoiling particle. The bottom panel shows the differential cross section where this m2 dependence cancels over 
most of the range. The insert shows any residual differences as a ratio of da/dQ for 2 different values of m2. 



respectively. Note that the scattering cross sections decreases for finite and small mass of the recoiling 
particle. This, in turn, leads to larger mean free paths A = 1/crp in QCD matter and smaller radiative 
energy loss. This is in contrast to the result of Ref. [43]. The reason for this difference is that in [43] 
the general term "dynamical medium" was inaccurately used to describe a specific hard thermal loop 
approximation and the reported increase arises from the lack of magnetic screening. Lattice QCD results 
and non-perturbative arguments, however, suggest that magnetic screening effects may be present already 
at 0{g^T). Our results for a general finite-range interaction mediated by a massive vector particle allow to 
precisely quantify the effect of the medium recoil. For example, we found that more than 90% of the cross 
section comes from configurations where the jet is not deflected more than 15% from its original direction 
of propagation. 

A small difference in the magnitude of the scattering cross section exists between the full calculation 
and the often used analytic approximation a — 87va'^/9/i'^ to quark-quark scattering. We can re-express 



the exact differential cross section as a function of the transverse momentum transfer as fohows: 

da da ^ C2{R)C2{T) |i;(q^; mi, ms)^ 

dn ^ d2q^ (2^)2 • ^ 

In Eq. (3.4) C2{R), C2{T) are the quadratic Casimirs for the incident and target parton representations, 
respectively, and rf^ = 8 is the dimension of the adjoint representation. The above expression also defines 
i;(q; mi, 7712), which now depends on the jet energy and the masses of the scattering particles. It is 
easy to check that such a definition reduces to i;(q; 0, 00) = ^i^^2 , consistent with similar definition 
in [40]. The only subtlety is that one allowed value of q^ generally corresponds to two values of cos(0). 
However, the region — 1 < cos(0) < contributes ^ 0.1% to the cross section and we simply ignore the 
second solution. 

Finally we discuss which momentum region for the exchange gluon gives the dominant contribution to 
the cross section. The momentum transfer in terms of the final and initial jet momentum equals q = ps —pi- 
Writing this in laboratory frame in terms of the light-cone components, we get: 

g+ = Es{l + cos(0)) - 2Ei , q- = Es{l - cos(0)) , \q^\ = Es sm{0) . (3.5) 

In the formula above we assumed that mi = and 1712 is arbitrary. As we can clearly see from Figure 1 
the cross section is dominated in the forward direction. Thus we assign to the leading region for the cross 
section the following power counting: 6 ^ X, £"3, Ei ^ 1. We immediately can see from Eq. (3.5) that this 
power counting corresponds to the following leading momentum region, depending on the mass 1712^'. 

q [A^, A^, A], if 1712 ^ A^, i.e. static source, (3.6) 
q ^ [X^, A^, A], if 1712 ^ A"*^, i.e. soft source. (3.7) 

To summarize, we investigated in detail the kinematics of jet-medium interactions to determine the 
feasibility of an effective theory where these interactions are mediated by Glauber gluons. We showed 
that for the static source, which we use in sections 5-9 below, the cross section is indeed dominated by 
the Glauber momentum region. We also calculated the exact dependence of this scattering cross section 
on the energy and mass of the incident and recoiling particles. This cross section will be implemented 
to go beyond the static scattering center approximation. We note, however, that this approximation is 
remarkably good and within 15% of the exact result even for low energy jets. 

4. An effective theory for jet propagation in QCD matter: SCETq 

As we have seen from the previous section, the effective theory of jet interactions in matter has to contain 
the Glauber mode, which carries the exchange momentum between the incident parton and the QCD 
medium with scaling q ^ [A^, A^, A]. Such mode is absent in SCET and we have to modify the theory by 
including it. One possibility was considered in [35] and later on used in [36] to study the multiple collisional 
interactions of jets. There are a few differences between our approach and these references. First of all, 
we write down the Glauber term directly in the momentum space as an effective potential [44], similar to 
the non-relativistic QCD (NRQCD) potential term [45]. Secondly, we consider a static source of Glauber 



^Such sensitivity of q'^ on the recoil mass arises from the fact that E^—Ei — 7712 — y^m^ + qr^, which is energy conservation 
equation. Thus, for heavy mass 7712 1 the recoil energy is negligible E3 — Ei ^ and for the lighter mass 7712 ^ A it is 
comparable to the transverse momentum transfer: Es — Ei ^ X. 



gluons, whereas in [35] the source was a massless cohinear field. We are motivated by the physical picture of 
nucleons or massive quasi-particles as sources of these Glauber gluons in nuclei and non-Abelian plasmas, 
respectively. We work in the rest frame of nuclear matter and also include the collinear gluons into the 
interaction Lagrangian with Glauber gluons. 

In a different context, namely the Drell-Yan process, it has been shown that the Glauber mode has to 
be added to SCET for the consistency of the effective theory [46]. Having formulated a consistent effective 
theory SCETq , it would be interesting to revisit the Drell-Yan factorization in this effective theory and 
understand the cancellation of Glauber gluons in inclusive Drell-Yan cross section from the effective theory 
point of view. In traditional QCD this cancellation was derived in Refs. [47, 48, 49], however it has not 
been addressed yet in effective theory methods to factorization. 

4.1 The SCETg Lagrangian for different sources and gauges 

Consider a quark or a gluon propagating in the positive light-cone direction n in QCD matter. In this 
subsection we will derive the effective Lagrangian of SCETq , which describes the interaction of our prop- 
agating jet with the source of the Glauber gluons. We consider three types of sources of Glauber gluons. 
The first one is a collinear field propagating in the n direction (considered in Refs [35], [36]). The second 
one is a (initially) static nucleon/nucleus or a massive quasi-particle. Interestingly, this source can be 
adequately described as a heavy quark effective theory (HQET) current. Finally, the third type of source 
that we consider is a soft parton of momentum p ^ [A, A, A]. For each of these three sources we consider 
three gauges: the gauge and two distinct light-cone gauges A'^ = and A~ = 0. In order to verify 
the symmetry properties of the Glauber interaction term we include the source fields into our Lagrangian. 
However, for practical calculations we integrate out the source fields as well as the Glauber gluons and 
present the resulting Feynman rules. 

We start from the vector potential as a function of the QCD current of the source and the gluon 
propagator. Our method to derive the scaling for the vector potential Aq(x) created by the Glauber 
gluons is same as in [35]: 



where we have written down only the part of the vector potential that contributes to particle-particle 
scattering, and the three remaining combinations involving anti-particles are omitted. In each of three 
gauges under consideration, the gluon propagator is equal to: 




(4.1) 



Expanding the propagator and the fermion field in the momentum space we get: 
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0) gauge 
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T- ^J'■ 
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Soft source 
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Ai/2 
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TqqAc 


[A^A^A^] 
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^1 

(n ^ n) 


[A2,A2,A3] 

pM 
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^1 


[A, A, A] 

^1 
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A^'{x) 

TqqAo 


^2 

(n o n) 


[0,A2,A] 

vn/iz^A 


[0,A,1] 

vn/iz^A 
^2 

^ 4 


A- = 


A^'{x) 

TqqAo 
FggAc 


[A2,0,A] 

pM 

^ 2 

^3 

(n ^ n) 


[A2,0,A] 

pM 

^ 2 

v^/ii/A 
^3 

pM 

^ 3 


[A, 0,1] 
p/i 

^ 2 

v^/iz/A 
^3 



Table 1: Summary of the scaling behavior of the Glauber gluon source ingredients, the Glauber vector potential and 
the Feynman rules for the newly constructed theory SCETq in the covariant gauge and two different light-cone 
gauges = 0, A" = 0. 

where the arrow in each line indicates that the term proportional to vanishes in Eq. (4.2) because the 
external source particles are on-shell: q^u(p-^ q)^yu{p)—{). In Eqs. (4.3) - (4.5) we also use a finite range 
interaction mediated by a vector field of mass ji. 

Using Eq. (4.2) and Eqs. (4.3) - (4.5) we obtain the scaling formula for Aq{x) for all cases. The 
scaling of the creation operators is derived from the anti-commutation relations in momentum space by 
allowing the external momentum of the source particles in each of three cases under consideration to be 
correspondingly: collinear, soft^ and soft. We also note that the vector p is collinear, while q has Glauber 
scaling for collinear and static sources, and q has scaling given by Eq. (3.7) for the soft source^. As a result 
we complete in table 4.1 all entries for the scaling of Aq{x). 

Having determined the scaling of the vector potential created by the Glauber field for the cases of 
interest, we now derive the effective theory Feynman rules for the coupling of the energetic jet to the 
Glauber field. We read off the Feynman rules from the usual SCET Lagrangian Eqs. (2.1) - (2.4) by 
treating the vector potential created by Glauber gluons in the covariant derivative as a background field. 
All Feynman graphs between collinear and Glauber gluons contained in Eq. (2.2)^ can be found using 
derived Feynman rules in the background field method from Ref. [50]. In addition to these rules, one has 

^For the static source the external momentum equals p — mv + /c, however pi — P2 — ki — equals the difference of two 
soft momenta. 

^Throughout this section we refer to both scalings in Eq. (3.6) and Eq. (3.7) as Glauber gluons, however it should be clear 
from the context which one is used when. 

^ Those rules include couplings of a Glauber gluon to two and three collinear gluons, and two Glaubers with two collinear 
gluons. We present in the table only the first vertex FggAc • 



to apply the specific power counting of the vector potential, which we derived for each source and gauge. 
As a result, we fill in the table the vertices FqqAG^TggAG for each source and for each gauge. Finally, we 
include in table 4.1 the Feynman rules Fg for the interaction between the source fields and the Glauber 
gluons. This is achieved by noting that we can view the jet moving in the direction of n as a source of 
Glauber gluons for the target fields 77, 77. Thus, the Feynman rules F^ can be found by using our scaling 
rules for Aq{x) with the collinear source in the n— direction^. Also, note that for the collinear source we 
use SCET current in the n direction, for the static source we use the HQET current and, finally, for the soft 
source we use the unexpanded vertex 7^ consistent with the soft quark interaction with the background 
field. As a result we fill in all the elements of table 4.1, where we have defined F1...F4, XI1...XI3 as follows: 



^2 



^ 1 




^ 2 


n-p 2 ' 


^ 3 


= igT" v^" , 
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(4.6) 

(4.7) 

(4.8) 
(4.9) 

(4.10) 



(4.11) 



yr^fiuX^abc 
^3 



9f 
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+gf{p'i+pT) 



{p- -2p'-)+Pl-2pr]+gT [\ip'- 



2p- 



+ p'l 



2pi 



(4.12) 



The derived rules allow us to write down the effective Lagrangian of SCETq • As a result, we obtain 
the following interaction term between SCET collinear fields and the vector potential Aq{x) of the Glauber 
gluons: 



>CsCETg (Cn, An, Aq) = ^CsCExlCn, ^n) + (Cn, ^n, Aq) , (4.13) 

Co (U^An^Ac) = 5^e-(^-^> (wr^;;Aof - K^A^ Ky), (<p) J ^g.A^) • (4.14) 

Depending on the gauge and the source, the vertices and the vector potential are different and are provided 
in the table above. The Lagrangian of this form for the collinear source in and A~ = 0^ gauges was 
derived in [35] and agrees with our expressions for corresponding two entries for FqqAc in table 4.1. We 

^Note that in the table we have derived the Feynman rules generated by the collinear source moving in the ft direction. 
However, since our target is a collinear current in the n— direction, its effect on the source can be derived from our table by 
reversing the n ^ n in the collinear source column. 

^In order to avoid confusion we note that in [35] the source was in the n direction while the target jet in the ft direction, 
thus our formulas agree with that reference if n ^ n as expected. For example in [35] the light-cone gauge A'^ = was 
considered, while it is analogous to our A~ = gauge. 



also note that for the covariant gauge and ^ = 1 our Feynman rule for YggAQi^O ~ disagrees with 

that of [36]. The corresponding Feynman rule from Ref. [36] contains only the first term in Eq. (4.10). 
However, note that with such a Feynman rule, the relation Eq. (7.4) below would be violated, which would 
lead to different results for single Born radiative loss (see section 6) calculated in the covariant and the 
hybrid gauges^. With our Feynman rule this inconsistency does not happen. 

Finally, in order to analyze the invariance of our Glauber exchange terms for quarks and gluons under 
the gauge symmetries of SCET, we rewrite Eq. (4.14) including the source fields (see Eq. (4.1)): 



where all the vertices for the target and the source are provided conveniently in table 4.1. In order to 
make this Lagrangian collinear gauge invariant one needs to dress the quarks and gluons with collinear 
Wilson lines Wn{x), defined in Eq. (2.12). As a result the Lagrangian that includes the Wilson lines can 
be obtained as follows: 



where WA^n{= Xn), t3n{A n) are the dressed collinear gauge invariant quark and gluon fields, correspond- 
ingly. In the next subsection we will show that Lagrangian in Eq. (4.16) is invariant under the soft and 
collinear gauge transformations of SCET. 

The derived Lagrangian of SCETq in Eq. (4.14) and Eq. (4.15) contains only interaction between 
a single collinear quark or gluon with a single Glauber gluon. However there are additional interactions 
between the collinear particles and Glauber gluons. For example, in the light-cone gauge the first term 
of Eq. (2.1) contains an interaction where the two Glauber gluons interact at the same point with a 
collinear quark line. The same Lagrangian in the same gauge contains the interaction where at the same 
point there are a collinear quark, a collinear gluon and a Glauber gluon. While we omit these terms in 
this section for brevity, their derivation is straightforward and the corresponding Feynman rules are listed 
in the appendix C. 

While we derived the Feynman rules of SCETq for variety of sources and gauges, in the main part 
of this paper we will do calculations in the following cases. For the source, motivated by our study in 
section 3, which showed that at RHIC and LHC energies the recoil effect is negligible with accuracy better 
than 15%, we will use the (initially) static source. The interested reader will notice that since the physics 
results depend on the transverse momentum exchanges (Glauber gluons) between the projectile and the 
target and the jet-medium cross sections, they should not be sensitive to the components of the source in 
the n and n directions. 

As for the gauge choice, we consider three cases in this paper, which allows us to establish the gauge 
invariance of the broadening and the radiative energy loss results. The first two choices are the covariant 
and positive light-cone gauges. All relevant Feynman rules are derived in this section and summarized 
in appendix C. Note that in both cases the Glauber gluons and collinear gluons are quantized with the 
same gauge-fixing condition. However there is no problem in quantizing the collinear gluons in, say, the 
light-cone gauge, while the Glauber gluons are in the covariant gauge. This is our third choice and we call 
it the hybrid gauge. For practical purposes this turns out to be the most convenient gauge choice [40, 51]. 



Cg (Cn,^n,^) >Cg 



(4.16) 



This gauge choice is defined at the end of this section. 



Physically, such a choice is possible since the scattering part and the medium-induced splitting parts of 
the calculation factorize. From a formal point of view the hybrid gauge corresponds to the vector potential 
created by the source derived in the covariant gauge, and the SCET Lagrangian quantized in the light-cone 
gauge with the corresponding background field vector potential. We discuss in details the Feynman rules 
in this gauge in section 7.1 and in appendix C. 

4.2 Gauge invariance of Cq, 

In this subsection we show that the Glauber Lagrangian Cq is invariant under both collinear and soft gauge 
transformations in SCET. The collinear gauge symmetry is a simple consequence of dressing the collinear 
fields with the collinear Wilson lines. The fields transform under the collinear gauge transformation of 
SCET according to [11]: 



1. 



Uc Ul + -Uc iD^ Ul , 



9 



Xn 



(4.17) 
(4.18) 
(4.19) 



where W is the collinear Wilson line, which is defined in Eq. (2.12), and the square brackets in the last 
equation indicate that the derivative operator acts only within the brackets. Note that the massive fields 
hy do not transform under the collinear transformation, since ph ^ [1, 1, A] and Pc ^ [1, A^, A]. As we see 
from the definition of Cg{Xi ^5 ^^v)^ since it is explicitly built out of the gauge invariant collinear fields, the 
Lagrangian is invariant under collinear gauge transformation: 

^ . 10 7 \ collinear gauge transformation r» / ir, i \ /^^^\ 

Cg{x, hy) ^ > Cg{x. hy) . (4.20) 

Demonstrating the invariance of the Lagrangian of SCETq under the soft gauge transformation is 
slightly more involved. The soft transformation of SCET looks like [11]: 

in ^Vsin, Xn = WUn-^VsXn, (4.21) 



9 



(4.22) 
(4.23) 



where Vg = e*"''^^)^", such that dfj_Vs{x) = ©(A^). For shorthand notation, define Y = Vg. The quark part 
of the Lagrangian of SCETq then transforms into: 



soft gauge transformation 



'-kl 



{Y^T^Y\, {Y^T^Y)u = (r^r^n^n,) ^ {YlfT^gYgi) 



2N 



:Sil6. 



1 

2N 



rpC 



Similarly, for the gluon part of the Co in Eq. (4.16) we have: 



soft gauge transformation 



-yaa' rpa' -ybb' rjii 



Y = A^"'"' A"'''' Y'' YT^'yK YT^'y'^ 



> 3^"" y A 



Y 



TIC 
u ■ 



rpa rj-ii 



At^'^-'A^'^'if"-' 



b'c' 



Y 



(4.24) 



As a result, we obtain that the Glauber term in SCETg is invariant under the soft gauge transformation: 

r ( t2 u \ soft gauge transformation r ( n u \ I A ^r.\ 

Cg{x,B,K) >Cg{x,B,K). (4.25) 



5. Jet broadening 



In this section we derive the modification to the transverse momentum distribution of jets from an elastic 
in-medium scattering to first order in opacity. We consider a quark jet or a gluon jet interacting with an 
initially static fermionic center for definitiveness. The necessary effective theory SCETq Feynman rules 
were derived in the section 4 and we use the covariant gauge for our calculation. The amplitudes that we 
consider for the quark and gluon cases have the following form, respectively: 

A^^^ = (J|rxn(xo)e^^'''^(^^'^°+^^'^^^^)|p), (5.1) 
^fe) = (J,,|rer(-o)e^^^^K^--+^^^-^)|p), (5.2) 

where Xn and Bn are the gauge invariant quark and Gluon fields in SCET, and J is the underlying hard 
process that creates the quark or gluon jet. The C^^^ term simply generates the hard QCD process in 
question, which we take into account by effective Feynman rule (J| Xn(^o) |p) = Xn,piJ{p) ^^^^^ - In the 
next two subsections we calculate these amplitudes using the Lagrangian of SCETq and combining the 
single Glauber gluon exchange diagram with the contact limit of the two Glauber gluon exchange diagrams. 
The corresponding amplitudes are called single and double Born diagrams in the literature [37, 38]. Each 
interaction in medium can be considered located at a certain point x, and an integral over the Glauber 
gluon momentum is introduced. To keep formulas compact, it is convenient to use the following shorthand 
notation: 

d^i = ^ e"^^'^^ v{qi) , (5.3) 

where i is the index of the corresponding scattering center, v{q) = 27T6{q^)v{qj_), ^{q±) is defined by 
Eq. (3.4) and Sxi = Xi — xq. The delta function in v{q) makes the dq^ integral trivial, while the q~ integral 
still needs to be evaluated. The following identity arises once one works out the d'^q in the light-cone 
coordinates: 

d$i = d*^x^e^^"^^^ (5.5) 

where 6zi = 6x^. To also shorten the color notation we substitute the color generator by simply a and 
specifying the representation it belongs to. For example: 

(a)fl(a)r, =r«(i?)T''(i), (5.6) 
[a,bU = ir''^T%R). (5.7) 

We use R to denote the color representation of the incident energetic quark or gluon, and we use T to 
specify the color representation of the scattering center itself. For example, = 3 since we consider a 
fermionic scattering center, while dji — dp — for a quark jet and dn — dA — ^ for a gluon jet. 

For the quark (or gluon) propagator with the collinear plus Glauber momentum p — q, where p has 
scaling [1, A^, A] and q is the Glauber gluon with the momentum scaling [A^, A^, A], we define: 




Figure 2: Three lowest order diagrams contributing to a collisional in-medium quark interaction. Top: tree level, 
bottom left: single Born diagram, and bottom right: two Glauber gluons exchange diagram. The notation for the 
scattering centers is the following: 01 = [xi^qi^ (^i)i],(02 = [^2,<?2, (^2)^]- The contact limit of the amplitude A'^^^ 
is also called double Born amplitude. 



It is convenient to define the following quantity, since it directly appears in the expression in Eq. (5.8): 

Note that the quantity in the Eq. (5.9) depends only on qj^ and p and always leads to a q~ pole in the 
propagator Eq. (5.8) which is in the positive imaginary plane. 

5.1 Quark jet 

We start from the first diagram in figure 2, which represents the lowest order diagram for the matrix 
element in Eq. (5.1): 

4°^' = Xn,piJb)e^^^°. (5.10) 

For a k Glauber gluon exchange diagram one has a factor of ^ for each quark propagator and a factor 
of I for each vertex. These factors can be commuted through to the left and simplified: 

ItV-y (lf]-Y (511) 

which means that we can (and will) ignore this factors in the Feynman rules. Squaring the tree level result 
in Eq. (5.10) we get: 

E l4°^'l' = Tr(^J(p)J(p))n-p. (5.12) 

spin, color 

It is instructive to examine here the Lorentz structure of the part of the matrix element squared which is 
not explicitly written in Eq. (5.12). It can be represented completely generally as follows: 

16 

J(p)J(p) = 5^a,(2.)r . (5.13) 

i=l 

Here, F^ is a basis set of 16 matrices for the Dirac algebra, for example: {1,7^,7^,7^7^,0-^^ = |[7^5 7^]}- 
Note that a^-^ = e^-^^S^. By direct inspection one sees that only the a+(p)7+ part contributes to the spin 
averaged cross section. 



For a k Glauber gluon exchange amplitude we have^: 



A 



(k)q 
coll 



Xn,p 



m=l \ 1=1 J 



(5.14) 



where S^*^)"? is the non-trivial part of the amplitude given by the SCETg Lagrangian. With the notation 
defined above and using the Feynman rules derived in section 4 and in appendix C we get in the i?j gauge: 



B 



(0)g 



(5.15) 
(5.16) 
(5.17) 



We have absorbed all the factors v{qj_) and the phases into the conveniently defined differentials in Eq. (5.3). 
Next we need to perform the longitudinal integrals dq^, dq~ in the one and two Glauber gluons exchange 
diagrams and reduce the integration to only transverse components of q^j_. For the single Born diagram 
this is done by using the result for the longitudinal integral /^^^ from appendix D.l. The result is: 



(5.18) 



Similarly, using the results from appendix D.l for the integrals and its contact limit we get the 

following results for the two Glauber gluons exchange diagram and its contact limit: 



J{u;i2Szi+uj2{^Z2-Szi)) 



(-1) 



/ 



(2c)q 



hbi)R(bib2)Ti 



In Eq. (5.18)-Eq. (5.20), the inverse formation times cui^2 are defined as follows: 

uji = n{p, qi) , UJ2 = 0,{p, q2) , L0\2 = ^{v, Ql + Q2) 



(5.19) 
(5.20) 

(5.21) 



Finally, combining together the three amplitudes, squaring them and identifying the contribution to first 
order in opacity, we get: 



^ Tr|4°)^ + 4^)^. 



dndr 



Here, we have omitted the Re (^A^^^^^ A^^''^ term, since it vanishes because of the color trace: Tr(T^(i?)) = 0. 

Squaring the A^^^^ term has been performed in Eq. (5.12) and represents the squared matrix element of 
the production of the quark jet from the underlying hard process J. The same exact overall factor, the 
differential jet distribution d^Gq^g/d^p^, appears in all other terms in Eq. (5.22). We will drop this factor 
for brevity but will keep in mind that the results derived below should be understood as operators acting 
on the unperturbed by the medium jet distribution. Thus, the first term becomes simply unity: 



1 



(0)<?|2 



+ |4'^^|2 + 2Re [A\ 



l(0)<? 



A 



(2c)q 



+ ... 



(5.22) 



dndr 



(5.23) 



In our notation B 



{k)q 



B stands for "Broadening" . 



Squaring the single Born amplitude A^^^^ we get: 
dn (It 

(5.24) 



11^ r 



This expression can be simplified further by turning the sum over the scattering centers into a continuous 
integral. This approximation is valid if ^ 

Y: e-^-^- - N^^^^^^^ . (5.25) 

i=l J ± ± 

Using this representation and integrating over q'-j_^ we obtain a particularly simple expression: 

In Eq. (5.26) e~^^'^^^ indicates the shift in the transverse momentum of the initial jet distribution. The 
contribution of the double Born contact term to the first order in opacity cross section correction becomes 
after a similar averaging procedure, which sets + = in Eq. (5.20): 

Note that in Eq. (5.27) there is no net transverse momentum transfer to the jet and no momentum shift. 
Combining the contributions up to first order in opacity we finally get: 

— Tr f \Af ^r + 1 4^)'' 1 2 + 2Re (4°^^)^ 



dn dr 



(5.28) 



where we expressed the answer through the elastic scattering, which in the lowest Born approximation 
equals to: 

daeijR.T) _ C2{R)C2{T) \v{q^)\^ 

d'q^ - dA {2ny ' ^'-'^^ 

5.2 Gluon jet 

We repeat a similar calculation as in the previous subsection for the gluon jet scattering off a fermionic 
scattering center. The first diagram in figure 3 equals to: 

The square of this matrix element describes the hard scattering production cross section and appears in 
every term and is dropped for brevity below. For k Glauber gluon exchange with medium amplitude we 
use the following notation, similarly to the quark case above in Eq. (5.14): 

4oli!;, = ^''(P) [ n <la6 ^J''" (p~i2 A • (5.31) 

m=l \ 1=1 } 




Figure 3: Three lowest order diagrams contributing to tiie gluon jet in-medium interaction. Top: tree level, bottom 
left: single Born diagram, and bottom right: two single Born exchanges diagram. The notation for the scattering 
centers is the following: (g)i = [xi, 91, (ci)^], (g)2 = [X2,q2,ic2)j]- 



With this notation, using the Feynman rules of SCETq in covariant gauge we get for amphtudes for 
k = 0,1,2: 

(5(l).)^.,a. ^ jabc, ^P.(p _ ^)MP^^ ^ (5.33) 
(^(2)5)p.,a6 ^ fbrbc, ^^^^^^ |.p:<.(^ _ ^^^^^ -^^^^ _ ^^^J'lP: fab.c, (^^)^^ 

X Ef^HP -qi-q2,p- q2) {-i)A{p - gi - 52)^^^) • (5.34) 

In the Equations Eq. (5.33)-Eq. (5.34), A(p)^^^^ is a standard gluon propagator in the covariant gauge, 

and the T^^^ {p^p') is the effective Feynman rule for covariant gauge for Glauber gluon coupling to collinear 
gluons. As one can see from our appendix C (also from section 4) it reads: 

J^l^'ip, p') = gi'" n-p + n'' {p' - p)\ - n" {p' - pf^ ^ (n V + • (5-35) 

To work out the contraction of vertices for single and double Born diagrams, the following identities are 
useful: 

e^(jp) tl^ip -qi,p) N^^,p^\p - qi) = £±-{p± - qi±) + £^ n-p, (5.36) 
e,{p) tl'^ip - q2,p) Niy,{p - q^) Sf^^ {p-qi- q2,p- ^2) N^^i\p - q, - q^) 

P£±-iP±-Qi±-Q2±)n^, + in-pfe-t , (5.37) 



where we defined N^^^^ to be the numerator of the covariant gauge gluon propagator: 

Nl.y{p)^g,.- ^{1-0- (5.38) 

Note that even though we work in the covariant gauge, we are free to choose any polarization vector for 
the on-shell final state gluon. We assume for the external scattered gluon as well as its source the physical 
polarizations: 

n-e{p) = 0, p-e{p) = 0, (5.39) 
n-J{p-q) = 0, {p-q).J{p-q) = 0, (5.40) 

Equation (5.39) was the only assumption in deriving Eqs. (5.36, 5.37) above. If in addition one assumes 
Eqs. (5.40), then the first terms can be dropped in each of the Eqs. (5.36, 5.37) and only the transverse 



term survives. Under these conditions and using results for the same longitudinal integral as in the quark 
case, which are summarized in the appendix D.l, we get for the single Born amphtude^^: 



T,, 



i± 



(5.41) 



which differs from the quark case in Eq. (5.18) only by a color and Lorentz structure. Similarly for the 
two Glauber gluons exchange diagram and its contact limit we obtain the following results: 



r'^^V''^^Hc2)r,(ci)T, (yf) / d^id^i 



J{U012^ Zl+U02{^ Z2-S Zl)) 
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j(2c)<;\ 



[c2)T.{ci)T.{gT) j d^id^i 



X - 



(-1) , (5.42) 
(5.43) 



The final result for the cross section to lowest order in opacity coincides with the one obtained in the last 
subsection for the quark case, with the substitution of the gluon-quark cross section instead of quark-quark 
elastic scattering cross section. Combining the contributions up to first order in opacity we finally get: 



-V Tr ( l^f + + 2Re (a^^^^^^'^A 

dndr \ ^ / 
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daei{R,T) 
d'q± 



(5.44) 



Here, 



dalf{R,T) _ C2{R)C2{T) \v{q^) 



dA 



(2^) 



(5.45) 



In this case the R stands for adjoint representation (gluon jet) and T stands for fundamental representation 
(fermionic static center). 



6. Medium-induced bremsstrahlung 

In this section we use the Feynman rules of SCETq to derive the probability for an energetic quark to emit 
a gluon, induced by the jet interactions in QCD matter. This is equivalent to evaluating the differential 
distribution of the number of emitted gluons. We first present this calculation in the vacuum using SCET 
and later in the medium using the new SCETq Lagrangian. In each case we consider the covariant gauge 
and the initially static source. We also focus on final-state (FS) radiation. In the literature, such a 
calculation is typically done in the soft (emitted) gluon approximation. However, in SCET and SCETq 
dynamics, the leading interaction describes the collinear gluon emission, which will allow us to go easily 
beyond the conventional uj <^ E limit. We will perform this new calculation below in section 9, while in 
this section we will focus on taking the soft gluon limit and comparing to the previously derived results 
for radiative energy loss in QCD matter. 

A.en. = (J|Tx.(xo)e^/^^^(^"'"+^^'''^)|p,k). (6.1) 

^"^Note that in our notation the covariant gauge gluon propagator of colUnear plus Glauber momentum can be written as 
A(p — q)^^^^ — ^(R^){P ~ Q)^g{PjQ) Using this identity along with Eqs. (5.36,5.37) leads to the same results for the 
broadening of gluon jet as for the quark jet up to the color and Dirac structure. 
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Figure 4: Diagrams in SCET describing splitting of a collinear quark into a collinear quark and a collinear gluon. 



To study gluon emission, we start from the matrix element, Eq. (6.1), where J is the underlying hard 
process that creates the quark jet, Xn is the gauge invariant quark field, and k are the momenta of the 
final state quark and of the emitted gluon, correspondingly. Since in this section we consider only the 
case of the initial quark jet, we omit the quark index in the amplitudes below for brevity. The matrix 
element in Eq. (6.1) gets contributions from 0, 1, 2, ... Glauber gluon exchanges between the collinear quark 
and/or gluon and the sources in the medium. The first three correspond to vacuum emission, single Born 
amplitude and two single Born exchanges amplitude, respectively, and are calculated in the subsections 
below. To simplify the notation we write the n— Glauber insertion amplitude in the following form^^: 

^^"^'" = 5Xn,p(^n/ R{qi....,qn)^''^^^''ij{k + p-Y,qi}j e'^'^+P^^ e^{k) . (6.2) 
6.1 Obtaining the Altarelli-Parisi splitting function in SCET 

A large process is accompanied by bremsstrahlung even in the absence of in-medium interactions. 
Knowledge of the corresponding amplitudes is also essential for the evaluation of the interference effects 
between the different sources or radiation for jet production in the QCD medium. 

Calculation of the vacuum diagrams in figure 4 leads to the Altarelli-Parisi splitting function for the 
q qg process. This calculation has been performed in Ref. [52] in the light-cone gauge. We perform the 
same calculation here in the covariant gauge. We also demonstrate how in the small x = /p'^ limit the 
relevant radiation piece can be identified at the amplitude level. 

The calculation of the relevant radiative matrix element can be performed using different set of fields 
in SCET (see [53]): gauge dependent fields or gauge independent ones Xn and Bn- In the former 

case the Feynman rules contain Wilson line emissions, from the gauge invariant quark field Xn — W^nCn- In 
the latter case these diagrams are absent, since we do calculations directly with Xn- However, the difference 
in the second case is that the SCET Lagrangian is modified and contains explicit collinear Wilson lines 
Wn which changes the Feynman rules. The equivalence of different formulations was shown in [53] on the 
example of the quark jet function at one loop. We choose the first case, i.e. use fields ^n.^n to do the 
calculation and include the Wilson line graphs, while keeping the SCET Lagrangian free from any collinear 
Wilson lines. In order to avoid confusion, we note that for the external quark spinor we use Xn as before 
for the broadening, which is simply a matter of notation. 

We first recall that for the jet production amplitude we have: 



(6.3) 



spin, color 



^^In our notation stands for "Radiation' 



One notices from Eq. (6.3) that the contribution to the cross section comes from the part J{p)J{p) oc 7"*". 
For the case of gluon bremsstrahlung, the first ampHtude in figure 4 reads; 



\ n-[p-\-k) n-p n-pn-[p + k) J [p ky 

Even though we work in the covariant gauge, we are going to use physical polarization vectors for the 
emitted real gluon. We choose the following gluon polarization vectors, which are the only possible solutions 
that satisfy the two conditions = 0, k-£{k) = 0: 

enk)=(o,^-^f^,Si^y i = l,2. (6.5) 

For such polarization vectors the second diagram in figure 4 vanishes, since it is proportional to n^, and 
also the last term in Eq. (6.4) vanishes once dotted with the polarization vector. By contracting the 
polarization vectors from Eq. (6.5) with the SCET amplitude in Eq. (6.4), we get-*^^: 



el , (6.6) 



A± = k±{l -x)- p^x , (6.7) 

where x is the fraction of energy taken by the emitted gluon k^ = xp^ , p^ = {1 — x)pq. To derive Eq. (6.6) 
we also assumed that the external quark and gluon are on-shell. Squaring this amplitude, averaging over 
the colors of the radiating quark we obtain: 



l^l^J^^^of = iT^{t--PJJ (<-)'/ (<-)]) , (6.8) 

7° •-) ' 7° {Rf^ - 4C^^ (1 - - + f ) ^ (I)color (I)D.ac • (6.9) 

The fact that the expression in Eq. (6.9) is a Dirac scalar means that the cross section of splitting factorizes 
into the cross section of jet production times the Altarelli-Parisi kernel, which has a nice probabilistic 
interpretation. Setting the light-cone direction along the initial quark before splitting, i.e. = — we 
reproduce the Altarelli-Parisi kernel for q ^ qg splitting: 



= (1 - x)Tt (^po+ J(o) j(o)) X ig^c'F (^ " ^ + y ) ^ = l^-^-?!' >^ • (6-10) 



We note that the 1 — x factor is associated with the reduction of the cross section to have a high momentum 
final-state quark when a gluon is emitted. From Eq. (6.10) we can identify the radiative correction and 
supplying the one body phase space, written as dk^(Pk±/{27v)^2k^ ^ we obtain: 



dx(Pk± TT^ X k\ 

One easily recognizes in Eq. (6.11) the usual Altarelli-Parisi splitting probability for quarks (up to x — > 
1 — x). Taking the small-x limit one can conveniently reproduce the boundary condition [40] associated 



^^Note that in Eq. (6.6) all indices are contravariant, while for example terms in Eq. (6.4) contain contractions between 
covariant and contravariant vectors, like = p^j±i^ = 7l- 
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Figure 5: Single Born diagrams contributing to the medium- induced gluon bremsstrahlung Eq. (6.2). The notation 
for the scattering centers is the fohowing: = [xi^qi^ (^i)i]- 



with hard jet production that is subsequently used in the reaction operator approach to parton energy 
loss [40, 51]. 

We finally point out that in the soft gluon limit x <^ 1 the amplitude itself, Eq. (6.6), reduces as 
follows: 

Ajg^gg « e^f-° Xn,pJ (-5^^^^) = Ajg X M5^<^ , (6.12) 

and, just like in quantum electrodynamics in the soft photon limit, a radiation matrix element can be 
identified at the amplitude level: 

Ml^^ ^ -g ^^%^^ TV (6.13) 
6.2 Single Born amplitudes in SCETq 

In this subsection we derive the single Born amplitudes in the fully covariant gauge for the scattering off 
the initially static source. In figure 5 we list all the relevant single Born diagrams. In this gauge two 
additional diagrams A/^^A^ appear when the collinear gluon appears from the Wilson line W"^ . Feynman 
rule for such collinear gluon vertex is well known: 

rr(A:) = ^r; 



where k is the outgoing gluon momentum from the Wilson line. Although we use the gauge in this 
subsection, just as in the vacuum case above, we use the physical gluon polarization vector for the emitted 
gluon Eq. (6.5). With the Feynman rules of SCETq , derived in section 4, and the notation of Eq. (6.2), 



we get the following expressions for the amplitudes in figure 5: 



R\"^ + (p + fc)2 + ^/ (^i)fi(^i)^^^^g(P + ^'gi)' (6-14) 

X 7v(5^)(fc - iA,(p + k, Qi) r^""^ (6i)t. Sr (fc - 91, fc), (6-16) 

77/ • rC 

= , (6.17) 
4^^^'^ = . (6.18) 

Note that in the cohinear gluon vertices in diagrams A\^A2 we omitted the last term proportional to 
because after contraction with the polarization vector of our choice this term vanishes, since — ^. For 
exactly the same reason diagram A/^ vanishes. However, the reason why we ignored n^^ in diagram A3 and 
why A^ vanishes, is slightly more involved. The point is that both A3 and A^ have a common factor given 
by Eq. (5.36) with p ^ k. Since from this identity it is obvious that n^^ times this combination vanishes, 
we are allowed to omit this term in A^. For the same reason ^5=0. 

In order to reduce the integral d^iR^^^ to the d^i± integral we use the identity in Eq. (5.5). Also, 
substituting Eq. (5.36) into the expression for R^"^ makes it obvious that the entire dependence on 
appears through the propagators Ag(p^q). Using the form of this propagator from Eq. (5.8) we define 
the relevant longitudinal integrals l[^\l2^\lo^\ We evaluate these integrals in appendix D.2. Thus, using 
Eq. (5.5) and the expressions for I^i from appendix D.2, we get after taking limit x <^ 1 in R^l 3: 

J « i-i) (abi)Ribi)T,'^^^^ J d^ixe-o"^^! , (6.19) 

J d^iR^^^^'^'e^ik) « i-i) {ha)R (5i)r, J [l - e-°^^i] , (6.20) 

/ d^iR^^'^'^'^e^ik) « i-i) [a, h]^ (5i)t, / ^ ^f.^ ~ ^'^t^^ e^^'^^ [e— ^^^^ - ll . (6.21) 

J J l«± - <ll±) L J 

The two inverse formation times loq and loi are defined according to: 

uo = — X , ^1 = T . [b.ZZ) 

xPq xPq 

6.3 Double Born amplitudes in SCETq 

In this section we calculate all diagrams in figure 6 in the R^ gauge for the initially static source. Expressions 
for all of these diagrams are obtained directly from the Feynman rules of SCETq in the covariant gauge, and 
are straightforward, though lengthy. However, a nice compact relations can be found for these diagrams by 
relating them to previously calculated amplitudes R^^\ R^^'^ and S^-*^)^, B^^)^, B^^)^, 5^^)^, in the notation 
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Figure 6: Two single Born exchange diagrams (double Born diagrams in the contact limit) contributing to matrix 
element in Eq. (6.2). The notation for the scattering centers is the following: ®\ = [xi,qi, (&i)i], ®2 = [x2, 92, {b2)j]- 



of Eq. (5.14), Eq. (5.31) and Eq. (6.2). The result is as follows: 



= i?i^)^'"(p, k, q2) Bi'^'ip + k-q2, qi) , (6.23) 

{R^^Y'^ = Bi'^'ip, q,) R^^^^^ip - q2, k, qr) , (6.24) 

- ^s°^"'(^' ^ - - '^^^ K^T"''' '^^^ ' (^-2^) 

= 5«^(p, ^'^'''"(P - 92, A:, gi) , (6.26) 

= 92) Bi'^'ip + k-q2, 9i) , (6.27) 

[R'i^Y'" = 4')^'"(p, fc, 92) + - 92, 9i) , (6.28) 

(^rig)"'" = . (6.29) 



Next, we integrate over the longitudinal momenta and reduce the integrals to transverse ones. This is 
analogous to the procedure for the single Born diagrams. The only poles in come from propagators of 
collinear momenta plus Glauber momenta (given by Eq. (5.8)). We summarize the corresponding integrals 
in appendix D.2. Since the amplitudes ^7^3,9 vanish for physical polarization of the radiated gluon, we do 
not even consider the corresponding longitudinal integrals. Finally, using our results for the corresponding 

(2) (2) 

longitudinal integrals /];,••• /g from appendix D.2 and taking the x <^ 1 approximation, we get the 



following expressions: 
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(6.30) 
(6.31) 

(6.32) 

(6.33) 
"1 (6.34) 

(6.35) 
(6.36) 

In order to understand the lowest opacity contribution to the induced bremsstrahlung, one needs to combine 
the single Born diagrams computed in the previous section with the contact double Born limit of the two 
single Born exchange diagrams. The contact limit of two single Born exchange longitudinal integrals 
is derived in appendix D.2. Using these results we get the following contact limits (or double Born 
amplitudes) : 
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where loq and coi are defined in Eq. (6.22) above and W2,<^i2 are equal to: 
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(6.38) 
(6.39) 



All results in this and the previous subsections, derived in the framework of SCETq , agree with the soft 
gluon approximation previously derived in the literature [40, 51]. For example, to first order in opacity 
and without explicitly showing the integral over the position of the scattering center we find: 



dN9{FS) 
dk+d^k I 



daei{R,T) 



kUk^-q±r 
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(6.40) 

Note that in the equation above the superscript "g" stands for the radiated gluons and should not be 
confused with Glauber gluons, denoted by "G". However, we can go beyond that and calculate the finite-x 
corrections to single and double Born diagrams, similarly to the full Altarelli-Parisi splitting kernel, and 
not just its soft gluon limit. In section 9 below we derive analytical formulas for these finite-x corrections 
to radiative energy loss at first order in opacity. 



7. Gauge invariance of the jet broadening and the medium-induced bremsstrahlung 
results 

In this section we demonstrate that the single and double Born amplitudes calculated in the previous two 
sections are gauge invariant. As it is known on the example of SCET, the gauge structure of effective 
theory is more rich than that of a full theory. This is a simple consequence of having multiple modes 
for the gauge field. In our calculation we deal with two types of gluons: collinear and Glauber. Thus, 
we can gauge fix these two modes completely independently without changing any physical result. Since 
Glauber mode is an off-shell mode, it is integrated out from the theory and is presented in the form of 
the potential term in Eq. (4.15). Thus, the only gauge freedom for Glauber gluons is the choice of the 
propagator /S.^y{q) in our effective potential, which in principle can be arbitrary. The collinear gluon field 
on the other hand is a truly propagating degree of freedom, with the corresponding kinetic term contained 
in the SCET Lagrangian. For each collinear gluon one could choose a certain gauge-fixing term. 

In the previous two sections we considered the fully covariant gauge, in the sense that both collinear 
gluons are quantized in the covariant gauge, and also for the Glauber Lagrangian we choose covariant gluon 
propagator /\^y[q)R^. Below we consider two alternative gauge choices and demonstrate there equivalence 
to the previous results. First, we consider a hybrid gauge where the collinear gluons are in the positive 
light-cone gauge and the Glauber potential is in the covariant gauge. Second, we choose both the collinear 
gluons and the Glauber potential term in the positive light-cone gauge A^^g — 0. 

The equivalence of all considered cases can be formulated in the following way. All the diagrams in 
consideration are some combination of elastic scattering amplitudes with a real gluon emission amplitude. 
Since each of these two processes is gauge-independent, the resulting amplitudes, when contracted with 
the physical external gluon polarization vectors, are the same even when each of the gauges are chosen 
independently. This is equivalent to the statement that different modes for the gauge field in the effective 
theory can be gauge-fixed independently. 

7.1 Hybrid gauge = 0, R^{Ag) 

We start from the hybrid gauge because it is simpler from a practical point of view. In this case we consider 
the light-cone gauge for the collinear gluons, while the Glauber propagator in the potential is taken to be in 
the covariant gauge. The Feynman rules for this gauge are contained in figure 11. Note that the collinear 
Wilson line Wn = 1 is absent in this case. For example, for the single Born radiative energy loss case we 
only have the first three diagrams compared to the covariant gauge case in figure 5, while the remaining 
two are absent. Since the collinear gluons are in the light-cone gauge, the transverse gauge link at infinity 
could in principle add new Feynman rules, where Glauber gluons arise from the transverse Wilson line 
Tn (see appendix E for a brief review of this subject). However, in this case it doesn't generate such a 
Feynman rule for the following reason. The diagram for which the Wilson line generates a Glauber 
gluon that couples to the source vanishes in this hybrid gauge because the propagator of Glauber gluon 
doted with the transverse gauge field and the source term vanishes: A^^g'^^Vfj^ — 0. However as we will see 
from the next subsection this is not the case in the fully positive light-cone gauge, where this interaction 
plays an important role. This is another reason why the hybrid gauge is so convenient. 

For brevity, we consider explicitly only single Born diagrams for the radiative energy loss in detail and 
then quote results for the remaining cases. We have only the first three diagrams in figure 5 for this case. 



Using the Feynman rules from the appendix C we obtain: 

(^1 )hyb = [^^ + p+ + k+ + ^ J + + ^^^^^ '^^^^ + ^' '^^^ ' ^^-^^ 

(^Dhyb - <^i)r ibi)T. i A,(p, gi)i(a)fi [n'^ + ^'^^^^^'^'^^ + iA^^p + fc, g^) , (7.2) 

X5fV^i«''H^'i)T.n-fc. (7.3) 

While first two amplitudes immediately match those in the covariant gauge given in Eq. (6.14), Eq. (6.15), 
in order to work out the third amplitude, the numerator NjZ^ of the positive light-cone gauge gluon 
propagator should be inserted. The following identity is straightforward to verify using the Feynman rules 
at hand: 

e,ik)g'[^n-kN(^y'^Hk - qi) = e,{k) Ef (fe - q^.k) N'£^\k - qi). (7.4) 
Thus, we arrive to the result that = (^^123) • remaining two diagrams in the figure 5 

hyb 

are zero in the covariant gauge^^ and are simply absent in the hybrid gauge. 

Similarly, proceeding with the remaining single and double Born amplitudes for the broadening and 
radiation we verify that the hybrid gauge results match the covariant gauge results at the amplitude level, 
as they should: 



/hyb V J ' V /hyb V J R^ ' 

V /hyb V Jr. V /hyb V Jr. 



(7.5) 



7.2 Positive light-cone gauge = 



We consider the positive light-cone gauge in this subsection. The collinear gluons are treated in the light- 
cone gauge as well as the Glauber gluons. Similarly to the previous subsection, we consider in details the 
single Born diagrams for radiative energy loss and later on quote the result in all other cases. Clearly, 
the collinear Wilson line vanishes again in this gauge Wn — 1. However, the transverse gauge link (see 
appendix E) gives a new non-trivial Feynman diagram. In figure 7 we show the two new diagrams that 
appear in our gauge choice. The first one directly follows from the collinear part of the Lagrangian of 
SCET, once the background field with the light-cone vector potential scaling is added. It is summarized in 
figure 11 of appendix C. The second diagram in figure 7 arises from the Wilson line emitting a Glauber 
gluon which interacts with the source. The corresponding Feynman rule is derived in the appendix E and is 
summarized in figure 14. As one can see from this Feynman rule, it depends on the light-cone prescription. 
So does the light-cone (Glauber) gluon propagator. We will see in this section that these two dependences 



^^This is true only for the same choice of physical polarization vectors in the covariant, light-cone and hybrid gauge calcu- 
lations. If one chooses arbitrary polarization vectors, once squared and summed over polarizations the result is independent 
on them. However, in this case it would be impossible to verify the gauge invariance with the hybrid or positive light-cone 
gauge, at the amplitude level. 



Figure 7 : Two additional diagrams that appear in the hght-cone gauge for the single Born amplitude for radiative 
energy loss. Left: term arising from the SCETq Lagrangian in the light-cone gauge. Right: Wilson line 
contribution to the matrix element. The notation for the scattering centers is the following: 01 = \x\^q\^ (^i)i]- 



cancel non-trivially and the final answer in this case is identical to the covariant gauge calculation above 
for all light-cone prescriptions. Similar cancellation was found in [54] by introducing the T Wilson line to 
SCET and calculating the jet function at one-loop in the light-cone gauge. 

Using the Feynman rules of this gauge we evaluate the first three diagrams in figure 5, which are 
present in this gauge as well: 



p(l)/i,a 
^1,2,3 



^1,2,3 



+ ^^1,2,3 ' 



(7.6) 



where the correction terms Ai?, 



(l)/i,a 



AI 



sum up to the following expression: 

The two longitudinal integrals AI and A^+ are defined in the following way^^: 

dq^ ,„-A., 1 1 



(7.7) 



A/ = 



27r 



{p + k-qir [qt] ' 



A 



A+ 



J 27r 



[4] ■ 



(7.8) 
(7.9) 



The second integral is summarized in the table 7.2 and vanishes in the —ie light-cone prescription. However, 
the first integral is clearly non-zero in any of the prescriptions, since its value comes from the poles of the 
propagator denominator. Luckily the first diagram in figure 7 completely cancels this term: 



A/+(a6i)i^(6i)r,7±^i±A/ 



. (7.10) 



When the amplitude R^^"^ is added to the first three diagrams, the AI integral vanishes as expected, 
while the integrals stays. It vanishes only in the —ie light-cone prescription, exactly in which the 

transverse gauge link vanishes according to Feynman rule derived in appendix E. In all other prescriptions. 



^Note that q'^ = —q^ in this equations, set by the S{qi) in v{qi). 



Prescription 






+ie 


1 


1 


—is 


1 

k^—ie 





PV 


2 \^/c++i£ ' k+-i£ ) 


1 

2 


ML 


1 

fc++i£sign(/c~) 


1 
2 



Table 2: Prescription dependent integral appearing in the light-cone gauge. 



the contribution from the transverse Wilson line does not vanish and we include it by calculating the 
second amplitude in figure 7. We explicitly see that by including these Feynman diagram the prescription 
dependence cancels in all three remaining cases: PV, ML, which happens through pushing the q~ 

pole into the negative complex plane which leads to the zero integral. Thus, by including the T- Wilson line 
into the calculation we get the same result in light-cone gauge as in the covariant gauge for all light-cone 
prescriptions. To see this explicitly we apply Feynman rule for the single Glauber gluon exchange from 
quark T— Wilson line, using appendix E (see figure 14) and obtain: 



/ 



4") 



{ahi)R {h 



I 



71'' + 



J3+ + /C+ J3+ 



(7.11) 



where the longitudinal integral A^^ is defined according to: 
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f d<h 
J 2n 



Szl 
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qj + ie q{ 



(7.12) 



Using table E for the prescription dependence of C^^^ we verify that it vanishes in the —ie prescription-^^ 
and in all other prescriptions makes the At integral cancel exactly the contribution from integral 
above. Thus, in the presence of a transverse gauge link we verify the gauge invariance for all prescriptions: 



= . 



(7.13) 



Calculations in all other cases, like jet broadening and radiative energy loss in single and double Born 
amplitudes, are similar to case considered above. We quote the results: 



q,9 



q,9 



(7.14) 



Thus, for all light-cone prescriptions we get an unambiguous result, same as in the covariant gauge. If one 
ignores the transverse gauge link, one should employ one particular prescription, —ie in this case, in order 
to recover the correct result. For practical purposes this gauge choice is the most difficult from all three 
that we considered in this paper. 



^Note that in the —is prescription the integral Ay^+ vanishes as well. 



8. Identifying the in-medium jet evolution kernels 
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In the preceding sections we calculated the amplitudes for collisional and radiative jet interaction at a 
particular position xi and demonstrated the gauge invariance of the end results. It is clear that a direct 
calculation of all diagrams for any number of scattering positions x^^, xi^^ - • • is not possible. From the 
results at hand, however, we can deduce the effect of an in-medium interaction at the amplitude and cross 
section levels in momentum space and employ the resulting kernels to describe a number of collisional 
and radiative processes in cold and hot nuclear matter as solutions to algebraic recurrence relations with 
suitably chosen initial conditions [40, 37, 38, 51, 55, 56]. 

In the absence of long-range color correlations in the target, the relevant interactions at position to 
lowest order in as that can build up to a jet-medium cross sections (a total of two Glauber exchanges in 
the forward cut diagram) can be expresses as follows: 

(8.1) 
(8.2) 
(8.3) 

Here, a represents a set of relevant quantum numbers, such as longitudinal momentum, transverse momen- 
tum and color. The kinematic and color structure of the scattering is contained in the unit ( / ), direct {D) 
and virtual {V) operators, which evolve the amplitude Ai-^...i^_-^{a) of the propagating system. The unit 
operator indicates that there are no Glauber gluon exchanges between the projectile and the target at posi- 
tion Xi^ {in — 0). The direct operator indicates a single Glauber gluon exchange between the projectile and 
the target at position xi^ {in — 1). The virtual operator indicates two Glauber gluon exchanges between 
the projectile and the target at position xi^ {in — 2), i.e. in the contact limit. The set of amplitude indices 
in Eqs. (8.1)-(8.3), thus, encodes the complete history of jet interactions in QCD matter. Repeating the 
basic operator steps in Eqs. (8.1)-(8.3) any amplitude that includes parton scatterings inside the medium 
can be iteratively derived from the unperturbed jet production amplitude 

n 

Ai,...M = J] [5o,i„ + Kirr^bm + S2,iJ/m] Ma) . (8.4) 
m=l 

Time ordering is implicit in the above formula. The amplitudes A^^'"^"^ {p, c) are the complementary am- 
plitudes to Ai^...i^{a) given by 

n 

A'^-'-{a) ^ 4{a) n [So,^yi + Si,i^bl + S2,^J\ . (8.5) 

m=l 

The differential jet or radiative gluon distribution, depending on the problem, can be expresses as a sum 
over the interactions in the medium dN{a)/dPS — '^^=QdN^^\a)/dPS. Each contribution with a fixed 
number of interactions n can, in turn be represented as 

2 2 

dN^''\a)/dPS (X A''•••'-{a)A^,...^So^) oc Tr • • • A''-'-{a)A,,...M . (8.6) 

The trace is over any uncontracted color and spin/polarization indices. Using Eqs. (8.4), (8.5) we obtain 
a simple recursion identity which relates dN^'^^ / dPS to dN^'^~^^ /dPS through the reaction operator R 

dN^^^/dPS oc A''-'-'RnAi,...i„_, , Rn = blbn + K + t>J . (8.7) 



Thus, the most important step in investigating the effects of the medium on jet propagation is the identifi- 
cation of the reaction operator R. We note that in problems that involve long-range coherence effects, such 
as radiative gluon re-interactions in the QCD matter, the recurrence relations may be inhomogeneous. 



oc oo oo 




Figure 8: Three contributions to the cross section, associated with the application of the reaction operator, 
correspond to the three t = oo cuts. 



8.1 Reaction operator for coUisional interactions 

With the direct and virtual interactions of a jet calculated in section 5, we can identify the form of the 
reaction operator for collisional interactions acting on a jet of momentum p 



Rn 



dZnp{Zn) / 



Zn-l 



dcJel 



cre/5^(q±n) 



Here, L is the thickness of the tareret, h — i Vn , is the 2-dimensional impact parameter operator coniuerate 
to the t._ n,„„«ntun, p »!t,ng to the "ght and b* = .ts He.m.tian eonj.gate aoting to 

the left. In this paper we explicitly proved the gauge invariance of this reaction operator. Our results are 
not surprising since R is expressed in terms of physical quantities, such as cross sections and densities and 
kinematic modification to the initial jet distribution. 

The contribution to the jet transverse momentum distribution from n medium interactions can be 
written as: 



d7VW(p^) 



n 



dzj 
A 



^ d(Jel{Zi 



(ppj 



(8.9) 



where A is the mean free path of the jet in the medium. Summing over all fixed-n contributions that 
Eq. (8.9) specifies, we obtain the most general result for the broadening of jets that propagate and interact 
in strongly-interacting matter. 

Simplifications of the sum over the contributions given by Eq. (8.9) can be obtained for special cases. 
For example, for a uniform density the integrals along the path of the jet can be performed immediately. 
With the notation x = ^/A, the result reads: 
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This is the partonic version of the Glauber multiple collision series where the jet interaction in the medium 
are described by a Poisson distribution and the momentum distribution is modified by the normalized 
differential scattering cross section [37, 38] . 

A final simplification can be achieved if one recognizes that the individual jet-medium interactions are 
of finite range r = such as the one discussed in sections 3 and 4. Let us take the initial jet distribution 
to be in the positive light-cone direction, (iA^(p±)/d^p± = 5^(p±)- The Gaussian approximation can be 
best understood in impact parameter space where the Fourier transform of the normalized scattering cross 
section reads 



/ ^ J (27r)2 7r(q2_ + ^2)2 4^2 up ; \^ 2 ^7 ^ ^ 



Here h — |b| and in the quadratic term in Eq. (8.11) the log 2/(1.08 6) multiplicative factor has been 
absorbed into an approximately 6-independent constant ^. Fourier transforming back to momentum space 
we obtain: 

d^P± J (27r)2 27r XM^C 

The resulting distribution is of two dimensional Gaussian form and has a width of 2x1^^^^,, i.e. (pj^) = 
2XM^C- It should be noted that the Gaussian approximation Eq. (8.12) is only applicable for small 
transverse momenta. It misses the Rutherford scattering power law 1/p^ behavior for p^ > 2xM^- The 
general solution the problem of jet multiple scattering is given by Eq. (8.9) [37]. 

8.2 Reaction operator for medium-induced gluon bremsstrahlung 

The derivation of the reaction operator for radiative processes is significantly more complicated in com- 
parison to the one for collisional interactions. The first step is to reorganize the amplitudes obtained in 
section 6 and to identify the contributions that can be interpreted as an interaction of the parent parton, an 
interaction of the already radiated gluon (Sn) and, finally, a genuinely new source of radiation associated 
with this particular scattering (5^). 

We introduce the following notation for the radiative gluon inverse formation times: 

and transverse momentum propagators: 

XT _ p -q±ii -q±i2 q±im 

ki ' ^in^.■■■^r.) _ q^^^ _ q^^^ ' 

Further, we recall that each double Born exchange yields a factor —1/2 and the numbers of these exchanges 
in the amplitude and its complementary are: 

n—l n—1 

= NMi-i^-i) = E ^^^irn , = N,{A'^-'-^) = E ^0.^- • (8.15) 

m=l m=l 

Finally, the color matrices associated with the collisional interactions of the jet are denoted by 

TeKAi...i„_J = {an-iT-' ■ ■ ■ [aiY' , tUa'^-'-^) = [a^f-'^ ■ ■ ■ {an-if-'-' . (8.16) 



We first examine the single Glauber exchange and the direct operator reads: 



^.17) 



We note that the common factor to all diagrams 2ige • (• • • where the transverse polarization vector e 
contracts with the 2D propagators in Eq. (8.14), is not explicitly shown for brevity. It is understood in 
Eq. (8.17) that if^^^^ rotates the color of the radiated gluon d: if^^'^^d = [c, a^]. Next, we identify the 
virtual operator: 



Vn —\{Ca + Cr) — anSn — CLjiBn —CLnDn — \ {C A — Cr) 

evy^yj ^^ fi v-^n A J^CCLnCL / \ \ 
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Substituting Eqs. (8.17) and (8.18) in the definition of the reaction operator we find: 

Rn = {Dn - an)\bn " «„) " = (5, + Bn)\Sn + K) " Ca 
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^.19) 



The diagonal Bertsch-Gunion term contributes only for the first scattering n = 1 [40, 51]. The off-diagonal 
terms depend on the current I„, which in turn obeys a recurrence relation itself: 
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To summarize, the medium-induced bremsstrahlung depends sensitively on the boundary conditions - both 
at the amplitude and cross sections levels. 

Finally, we give one example for the complete solution to the final-state medium-induced bremsstrahlung 
for a jet produced in a large process. The boundary condition is represented by the amplitude: 
Jo = — 2z^ g^^o2;o ^ ^ associated with the real hard bremsstrahlung for x ^ 1. Rewriting the position 
integrals as integrals over the separation between the scattering centers Az^ — Zn — Zn-i and including the 
integrals over the momentum distribution of the jet-medium scattering cross section we obtain [40]: 
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To use unified notation above, we have to specify ^2 = B(^^i = B^. In the case of final-state 
interactions, ^ is the point of the initial hard scatter and zl = L is the extent of the medium. The 
path ordering of the interaction points, zl > ^j+i > zj > zq^ leads to the constraint Y17=i — ^L- One 



implementation of this condition would be Azi G [O^zl — X^j^i Azj ] and it is implicit in Eq. (8.21). The 
complete solution to the problem of medium induced bremsstrahlung in the x ^ 1 limit for three different 
boundary conditions can be found in Ref. [51]. With the proof of the gauge invariance of the reaction 
operator for collisional and radiative processes, we have also proven the gauge invariance of these results. 
As we pointed before, when the end result is expresses in terms of physics quantities, such as scattering 
cross sections, mean free paths and formation times, gauge-invariance can easily be recognized. 



9. Bremsstrahlung: beyond the soft gluon approximation 



In this section we calculate the corrections to the medium induced bremsstrahlung at finite splitting fraction 
X ^ 1. Effective theory Feynman rules allow us to do so easily. The advantage of the effective theory 
approach is that we have the appropriate interactions in medium at the Lagrangian level, which allows us 
to straightforwardly calculate any process of interest. Doing similar calculations, for example calculating 
radiative energy loss to first order in opacity keeping the full x dependance in the traditional approach 
would be more difficult, since it will require to do approximations at the level of Feynman diagrams. 



9.1 Incoherent radiation 

We first discuss the simpler case of gluon emission without the Landau-Pomeranchuk-Migdal destructive 
interference effects. Starting from expressions in Eq. (6.14) — Eq. (6.16) and without making approximation 
on X we identify three pieces of these matrix elements that are proportional to exp[iQi(z — zq)] (see 
appendix D.2 for details) as the ones that contribute to the Bertsch-Gunion amplitude - the QCD analog 
of the Bethe-Heitler radiation in electrodynamics. Squaring the sum of three diagrams and summing over 
the physical polarizations we obtain: 
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It easy to see that in Eq. (9.1) N/A± = dz p{z), where p is the density of scattering centers in the medium. 
Integrating over the path of the quark propagation we find that the differential spectrum of the incoherent 
Bertsch-Gunion bremsstrahlung can be written as: 
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Here, \g{z) = 1/ 
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p{z) and the initial jet direction (^^o — p + k)\\n. Very generally, one can 



express the radiative amplitude squared as follows: 
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In Eq. (9.3) 



A_L = (1 - a;)k_L - xpx = k_l |px=-k_L , 
B_L = (1 - x)k_L - x(p_L - q±) = k_L + xq_L |px=-k_L , 
C_L = (1 - a;)(k_L - q±) - xpx = k_L - (1 - a;)q_L |p_L=-kx ■ 
By taking the small a; <C 1 limit above and substituting in Eq. (9.3) we obtain the known result: 
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The general expression for the incoherent medium-induced gluon bremsstrahlung for a medium of constant 
density and length L without such approximations is; 
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Figure 9: Left panel: the differential gluon number distribution versus k± is shown for selected values of x = /p'^ 
for a quark jet that has undergone large scattering. A finite in-medium mass meff = 1 GeV regulates the collinear 
divergence. Right panel: the same spectrum for the incoherent medium-induced Bertsch-Gunion radiation. We have 
used L/Xg = 1 and for the x ^ also shown an analytic calculation in the massive scattering center limit for 
comparison to the exact calculation with TTimed. 

= 1 GeV. 



Finite-x effects are illustrated in figure 9. The left panel shows the Altarelli-Parisi real gluon differential 
spectrum Eq. (6.11) for selected values of x = 0.1, 0.2, 0.5, 0.7, 0.9 versus k±. We use x = 0.001 to 
simulate the small x-independent limit and represent the result by a solid line. The normalization is fixed 
by the choice of a quark jet and as — 0.3. A finite effective mass rngfr. = 1 GeV to simulate medium effects 
is implemented as A;^ ^ A:^ + m^^ and regulates any collinear divergence. The dashed lines show that the 
finite-x corrections can be as large as a factor of 2. The kinematic bound for the transverse momentum of 
the emitted gluon in large jet production k± max — ^x(l — x)Q'^ is not explicitly shown in the figure. 

The right panel of figure 9 shows the medium-induced differential gluon spectrum for the same values 
of X in the incoherent Bertsch-Gunion limit. We have chosen L/\g — 1 in Eq. (9.8) to facilitate direct 
comparison to the Altarelli-Parisi case. In our calculation we used the exact form of the 2^2 scattering 
cross section Eqs. (3.2), (3.3), (3.4), including the finite mass of the medium particles, the medium recoil 
and the kinematic bounds on max- Having chosen a quark of energy E — 100 GeV, we note that in the 
Bertsch-Gunion case large-x effects are even more pronounced than in the Altarelli-Parisi case. 

It should be noted that if an approximate form for the dependence of the normalized differential 
scattering cross section is employed: 

Gel d^q± 7r(g^ + M ) 

the remaining integral in Eq. (9.8) can be performed analytically. This form is motivated by the infinitely 
massive scattering center approximation. One uses a Feynman change of variables developed for the 
calculation of loop diagrams. We here quote the final result for the soft gluon approximation x <^ 1: 
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where Ai^2 are the solutions of the quadratic equation: 



Xmls + (1 - A) /i^ + A(l - A) ki = . 



(9.11) 



The analytical formula for the finite x case is also calculable, though slightly more involved. To illustrate 
the differences between the exact form of the 2^2 scattering cross section and the approximate form 
given by Eq. (9.9) we show the analytic small-x gluon distribution versus in the right panel of figure 9. 
We note that the differences between the exact and approximate calculations are very small and appear 
only at high A:^ - a part of phase space that does not contribute significantly to the medium-induced energy 
loss. 

9.2 Final-state radiation and the Landau-Pomeranchuk-Migdal effect 

In this subsection we calculate the combination of single and double Born amplitudes to first order in opacity 
by keeping the finite-x corrections. This evaluation proceeds analogously to the one for the incoherent 
Bertsch-Gunion limit. Directly from the Feynman rules of SCETg , derived in this paper, we find the 
combined squared amplitude from single and double Born diagrams equals to: 
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In the equation above F^, Gi are form-factors, which we give below, and q are the color factors: 



ci = C2{RfC{r)dR. 



C2 = C2{R)C{r) [ C2{R) - -Ca ] dn . 



We recall that in our example C2{R) — Cp^ (Ir — Nc and C{r) — |. The two operators that appear are 
the identity operator in the Dirac indices I and S^, which is the third component of the spin operator for 
the fermion: 

'^3 
as 
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The first term in Eq. (9.13) allows for the factorization of the medium-induced radiative corrections and 
the hard scattering cross section even in the large-x limit. The second term in principle does not vanish. 
If we, however, recall the decomposition of J{p)J{p) in the Dirac algebra basis, see Eq. (5.13), it is easy 
to verify that there must be a non-zero pseudo-vector contribution for this term not to vanish identically. 
Therefore, one might expect further corrections for processes with electroweak boson exchanges. For QCD, 
medium-induced radiative corrections always factorize. A simple exercise for the interested reader is to 
show this on the example of inclusive tree level jet production. 



The form- factors for the single Born diagrams in Eq. (9.13) are given by the fohowing expressions: 

(l-x+^y\/3,\' + \(3,\') , GP = 2(x-^)lm[(/3fr/3f + (/3fr^^], (9.15) 
F|B= (i-:, + ^^(/3^./3* + /3^./32) , Gf = 2(x-^) Im[(/3f)*/3f + (/3f)*/3|] . (9.16) 
On the other hand, the double Born diagrams contributions can be written as follows: 
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F^^^ (l-^+y)(ao-72 + «S-72), = 2 - Im [(7!)* + (ag)* 7^] ■ (9-18) 

In the expressions above the vector ao appears from the vacuum diagram, the vectors /3i 2 appear from 
the single Born diagrams and appear from the double Born diagrams calculation. They all are given 
in the equations below: 

2 A I 

c^o = ^, (9.19) 
/3i=(^-^)e^"^^^ + ^e^^3.., (9.20) 

/32 = - (^ - ^) e^"^^^ + ^ e^"^^^ - ^ e^"3.. , (9.21) 



7i = -^+(^-^) e^"^^^ + e^"^^^ , (9.22) 



3A, /2A, 2^\ ,^^5 2£i 



72 = -Tf- - -Tf- f- ' ^ e'''' , (9-23) 

(9.24) 



where the vectors A±, B±, C± are same as in the previous subsection, the phase factors ^1,2,3 are same as 
in appendix D.2. The remaining propagators and inverse formation times D±, 0,4, appear in the double 
Born diagrams and are equal to: 

Z)_L = A_L - g_L = ^± - Q± lpx=-kx , (9-25) 
42 42 _ r)2 

PqX[1 — X) PoX[l — X) 

First, we should check that our result agrees with previous calculation in the soft gluon approxima- 
tion [51, 40]. For this we need to expand all quantities to lowest order in x <^ 1: 

A^^k^, B^^k^, C^^k^-q^, D^^k^-q^, (9.27) 

^ CJQ 5 f^2 ^ , ^ CJo ~ '^l 5 ^4 ^ '^O 5 ^5 ^ '^O ~ '^l 5 (9.28) 

Gp^O, GP^^O. (9.29) 

Note that in this limit factorization is exact even if J{p)J{p) contains a pseudo- vector component. The final 
answer depends on the form-factors F^, which we calculate below in the soft gluon limit. It is convenient 



to rewrite the a,/3^ 21 7i 2 terms of standard definitions in the literature^^ [51, 40]: 
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The corresponding x — )■ limit is particularly simple: 



(9.31) 

/3i = 2Si e*'^°''^ + 2Ci e*('^o-'^i)'5^ , (3^ = 2Hi - 2Bi e'"""^^ - 2Ci ^^0-^^)^^ , (9.32) 
7i = -Sf/i + 2Bi e'^'o^^ + 2Ci e^^'^o-'^i)''^ , 72 = 21^1 - 2Bi e''^^^^ - 2Ci S'^^''^^'^^^ . (9.33) 

Useful relations between these vectors, that help in deriving the expression for the form-factors Fi are the 
following ones: 

/3i + /32 = 2Hi , 7i + 72 = -Hi . (9.34) 
From the definitions in Eq. (9.15)-Eq. (9.18) we get: 
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Finally, using these equations we combine the single and double Born form-factors into the sum: 
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Thus, in the soft gluon approximation we get: 

(/B ^ pDB)^^^^ ^ _ (_i6B,.Ci) (1 - cos(6.;ifc)) . 

Taking into account the phase space factors, the color factors and the final-state coherent medium-induced 
emission contribution above, we find: 
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in agreement with Eq. (70) of [51]. 

Beyond the soft gluon approximation, the full result for the coherent medium-induced bremsstrahlung 
reads: 
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^Note that Bi is distinct from B± and Ci is distinct from C±. 



We leave the discussion of this new result and phenomenological applications to new RHIC and LHC 
experimental data [29, 30, 31, 32, 33, 34] for future work. We note however that in Ref. [57] an evaluation 
of the medium-induced energy loss beyond the helicity amplitude approximation found ^ 18% reduction 
in the mean medium-induced energy loss. Our plan for the future is to investigate the large-x radiative 
correction reduction effects at the most differential level. 

10. Conclusions 

In summary, we constructed an effective theory SCETq for energetic quark and gluon p ^ [1,A^,A] 
propagation and interaction in dense QCD matter. This theory is well-suited to calculations both in the 
quark-gluon plasma [25, 26, 27, 28, 56] and in cold nuclear matter [38, 55, 56, 58]. 

To construct this theory, we examined the relevant t-channel parton scattering cross sections and 
demonstrated that they are dominated by forward scattering. The corresponding momentum exchanges 
q [A^,A^,A] are approximately transverse to the direction of jet propagation and mediated by Glauber 
gluons. We also demonstrated that a fully dynamical treatment of the scattering centers in QCD matter 
leads to a small reduction of the scattering cross section and smaller medium-induced effects in contrast to 
early speculations [43]. The SCETq Lagrangian was shown to be invariant under soft and collinear gauge 
transformations. We derived the Feynman rules for this new Lagrangian in the covariant and light-cone 
gauges. Also we provided a third choice, which we call the hybrid gauge, when the Glauber gluons are 
quantized in the covariant gauge, while the collinear fields are gauge-fixed with the light-cone gauge. This 
choice provides us with the simplest form of Feynman rules and is useful from practical point of view. 

The new effective theory was used to evaluate the broadening and medium-induced radiation for 
energetic quarks traversing a region of hot/dense QCD matter. SCETq was shown to recover exactly the 
known results for the transverse momentum diffusion of particles in the strongly-interacting medium [37, 
38]. It allows to study the Molliere multiple parton scattering beyond the limitations of the Gaussian 
approximation - namely, the inability to describe the Rutherford power-law tails of transverse momentum 
distributions. For the case of inelastic quark interactions we derived the fully differential medium-induces 
gluon bremsstrahlung spectrum for both the incoherent and Landau-Pomeranchuck-Migdal suppressed 
cases. In the soft gluon approximation we obtained the kernel for the reaction operator approach to 
medium-induced energy loss [40, 51]. We also evaluated the large x = k^/E^ corrections to the medium- 
induced bremsstrahlung spectrum. Gauge invariance of the jet broadening and radiative energy loss results 
was demonstrated explicitly for the first time by performing the calculations in the covariant i?^, light- 
cone and hybrid gauges. On the example of an energetic quark, we found that in QCD the process- 
dependent medium-induced radiative corrections factorize from the hard jet production cross section. This 
allows us to write down perturbative jet and leading particle observables in heavy ion collisions as a 
convolution of the corresponding observables in the more elementary nucleon-nucleon reactions and the 
medium-induced collisional and radiative corrections specific to the process under investigation. Our results 
put jet quenching phenomenology in heavy-ion collisions on more solid theoretical grounds. 

In the near future it will be a high priority for us to implement the improved theory of parton propaga- 
tion and energy loss in matter in perturbative QCD predictions for analysis of the energetic particle and jet 
quenching data from the heavy-ion experiments at RHIC and at the LHC. With the calculations at hand, 
we expect to be able to reliably combine the process-dependent medium-induced radiative corrections with 
next-to-leading order perturbative effects [59]. We plan to extend these calculations beyond applications 



to jet propagation in the quark-gluon plasma [25, 26, 28] and also improve the accuracy in the evaluation 
of cold nuclear matter effects [58]. 

The derived Lagrangian of SCETq can be used to revisit the factorization of the Drell-Yan process in 
the effective theory and include the spectator interactions into the analysis. The importance of Glauber 
(or Coulomb) gluons for the Drell-Yan factorization has been addressed, and their cancellation in the 
inclusive cross section has been proved in traditional QCD approach to factorization [47, 48, 49], while 
similar understanding in effective theory method is still missing [44, 46]. 
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A. Light-cone notation 

The momentum of any particle can be conveniently represented in light-cone coordinates. A light-cone 
vector /3 is defined by the condition = 0. The two vectors that specify the positive and negative light- 
cone directions are = (1, 0, 0, 1) and = (1, 0, 0, —1), respectively. An arbitrary four-vector p then can 
be expanded in the light-cone vectors basis: 

(A.l) 

where the +, — components are defined as follows: p'^ = n-p, p~ = n-p. The four- vector p^ thus can be 
written in the light-cone basis as [p^,J>~,p^], where we will use the square brackets to emphasize that 
the light-cone notation is being used. For example, the light-cone vectors have the following coordinates: 
n — [2,0,0], n — [0,2,0], corresponding to positive and negative light-cone directions. Thus, the plus 
component is in the n direction and the minus component is in the n direction. The degrees of freedom of 
SCET are collinear quarks and gluons, with the following momentum scaling pc — ^ A] and soft gluons 
with momentum ps — [A^, A^, A^]. The Glauber modes that we consider in this paper have momentum 
scaling pg — [A^, A^, A]. 



B. Kinematics in the laboratory frame 

In section 3 we estimated the recoil effect in 2 ^ 2 scattering of a projectile particle with mass mi on 
an originally at rest target particle with the mass m2. An important feature of the laboratory frame is 
that the final energy of the scattered particle is a non-trivial function of the laboratory angle 9. This 
function can be found from energy and momentum conservation laws. The solution is given by: 

, , Pi7^ zb (7^ — l) r sin6^cos0 

E^-rn^ '^ ^\ I • B.l 

7^ — (7^ — 1) cos^ 



In Eq. (B.l) we defined: 



PI = ^ , (B.2) 

PI 

72 (1 - pf) + cot2 e , (B.3) 



where /3i is the velocity of the projectile particle in the center-of-mass (CM) frame and /3 is the velocity 
of this center-of-mass. Note that the kinematics is very different depending on the masses. If mi > m2, 
then both solutions in Eq. (B.l) are physical, while for mi < m2 only the solution is physical, while 
the "—"-solution is not. In our case we can safely assume that mi < m2. 

The reason for the difference discussed above is that when the projectile particle is heavier then the 
target one, then in the center of mass frame the projectile particle's velocity is smaller than the center- 
of-mass velocity in the laboratory frame, i.e /3i < /3. This means that no matter in which direction the 
projectile particle goes in the CM frame, it will always move in the positive hemisphere of the laboratory 
frame with respect to the initial momentum of the projectile particle. Thus, each direction of motion of 
the scattered projectile particle in the laboratory frame corresponds to two distinct directions in the center 
of mass frame. That is the origin of two solutions from Eq. (B.l). Alternatively, for mi < m2, we have 
/3i > /3 and each direction of the projectile particle in laboratory frame comes from a distinct direction 
of this particle in the center of mass frame, leading to the corresponding unique energy in the laboratory 
frame. 
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Figure 10: Feynman rules of SCET for the interactions between collinear quarks and gluons. 



C. Feynman rules of SCETq 

In this appendix we review the Feynman rules of SCET and also present all the relevant Feynman rules of 
SCETg for different gauge choices. We consider the initially static source case. 

First, we start with the SCET rules. Using the Lagrangian of SCET [9] given in Eq. (2.1) one finds 
the Feynman rules of interaction between the collinear quarks and gluons given in figure 10. These rules 
are given in the covariant gauge. There are additional vertices in this gauge when one has more than 2 
collinear gluons and two collinear quarks at the same point, which we omit here. The same Feynman rules 
as in figure 10 are valid in the positive light-cone gauge A^, with the simplification that all the 
terms vanish. Also, all the mentioned omitted vertices vanish in the light-cone gauge. 

Our approach to deriving the effective Lagrangian of SCETq was to treat the effect from Glauber 
gluons on the target jet as one of an external background field, which the source creates. Having determined 
the scaling of the vector potential, we then read off the Glauber term from the SCET Lagrangian as it is 
given in Eq. (2.1) with a trivial addition of an external vector potential to the covariant derivative. We 
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Figure 11: Feynman rules of SCETq involving Glauber gluons. The first two vertices are given in the covariant 
gauge. Vertices from three through seven are for the light-cone gauge = 0. Last two vertices are for the Hybrid 
gauge. Note that the last three vertices in the light-cone gauge are power-suppressed in the covariant and Hybrid 
gauges. 



present in figure 11 the resulting Feynman rules. Note that we include the propagator of the Glauber gluon 
into the Feynman rule as well as the color and overall Dirac structure from the static fermionic source that 
we consider. 

In the main part of the paper we consider three gauge choices. First, we used the covariant gauge in 
the sections 5, 6. Second, we used the light-cone gauge and third was the hybrid gauge in section 7. In 
the first case we used covariant gauge for both pure SCET and the Glauber interactions given in figure 11. 
In the second case we did the same for the light-cone gauge. Finally, for the hybrid gauge, we used the 
light-cone gauge for SCET gauge fixing, while using the covariant gauge for Glauber gluons. In all three 
cases we found the same results for physical quantities, as expected. 

D. Longitudinal integrals 

In this appendix we define and calculate all longitudinal momentum integrals that are required for the 

(n) 

evaluation of jet broadening and radiative energy loss. Our notation is /m , where n indicates number 
of Glauber gluon exchanges, and m is the index corresponding to the Feynman diagram in question. We 
perform each integral exactly, without assuming a soft gluon approximation for radiative energy loss. 
However, from the exact results we can easily take the soft gluon limit. 

D.l Jet broadening 

The longitudinal momentum integral that appears in the single Born diagram for quark or gluon jet 
broadening equals to the following expression: 

^ f dq^^^q^6z, A{p,qi) = / ^e'^^^'^ ^ = -ie'^^^'^ , (D.l) 

where the inverse formation time uji is defined according to Eq. (5.8) and Eq. (5.9): 

The double Born longitudinal integral looks like: 

= (-z) J ^^'^^ (J((^i2-qo)Szi+q~Sz2) /^^[p^ g^) = _Qi{oJi25zi^uj2{5z2-5zi)) ^ (D.3) 

where inverse formation time (jJi2 equals: 

c.:. = + g.) = p- - ^P^-^^^- ^^^)'-'^ , (D.4) 

and (jj2 is identical to uji up to qi ^ q2- In deriving Eq. (D.l) we used the fact that ui is in the upper complex 
plane, which we should choose, since Szi > 0. Similarly in first and second step in deriving Eq. (D.3) we 
used that and (jJ2 are in the upper q~ complex plane. In the second step we also used that 6z2 > Szi, 



which is true for second order in opacity diagram (we use time ordered notation zq < zi < Z2 < ...). Next, 

(2) 

we calculate the contact limit of the integral /]; , when Sz2 = Szi + 0: 

= -ie^'-^i^^^i / ^ = ^ (ln(oo - UJ2) - ln(-oo - ^2)) = --e''^^^^'' . (D.5) 

J ZTV CJ2 — ^2 ^ 

As one can see, in the second step in deriving Eq. (D.5) we cannot use Cauchy's theorem to perform 
the remaining q2 integral, since the boundary term at infinity does not vanish. Instead, we perform the 
integral directly. We find: 

it^ ^ll[^HSz2^Szi) . (D.6) 

D.2 Radiative energy loss 

The longitudinal momentum integrals that appear in the first three single Born diagrams on figure 5 are:^^ 

/f^ = / ^ ^aiP + k, qi) , (D.7) 

^(1) ^ I ^g^.r^-i A,(p,gi)A,(p + A;,gi) , (D.8) 
^(1) ^ I ^ ^i,-s., ^^(^^ ^ _ 

These three integrals are functions of three inverse formation times, all given by the poles in Eq. (5.9) of 
propagator in Eq. (5.8): 

= nip + k, ^p- + k-- ^^^-^'^,7;/-" , (D.io) 

^2 = n{p, q,) =p-- - '^'f - , (D.ll) 

= nik, gi) = k-- - 91^^)' - . (D.12) 

The first integral in Eq. (D.7) we already derived in Eq. (D.l), where we need to substitute cji ^ Qi. 
The integrals l2^^ and I^^^ can be reduced to the same integral in Eq. (D.l) by using the following trivial 
identity, which follows directly from definitions in Eq. (5.8) and Eq. (5.9): 

Ag{pi, qi) Ag{p2, qi) = r ^rf — . (D.13) 

Note that in Eq. (D.13) the numerator depends on g^, while the denominator does not. Thus, integrals 
/2^^ and /g^^ become a combination of two integrals in Eq. (D.l). As a result, we get the following exact 
expressions: 

= -ie'^'^'^ , (D.14) 

= ' fe^^i^^ - e'^^^'A , (D.15) 
ill ~ ^^2 ^ ^ 



11'^ = _ [e'''^'' - e'''^''^ ) . (D.16) 



''The remaining two vanish when one takes the physical emitted gluon polarization vector. 



The soft gluon emission approximation corresponds to x <C 1 and is of interest in sections 5 and 6. In this 
hmit the inverse formation times f^i,2,3 reduce to only two non-trivial functions uq^cji. The corresponding 
approximate results read: 
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Next we move to the integrals in two Glauber exchange diagrams and their contact limits. The 
longitudinal momentum integrals that appear in the first six diagrams in figure 6 equal: 
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The corresponding integrals contain seven inverse formation times, which later after taking the contact 
limit and averaging over the medium states reduce to only two inverse formation times. We define the 
following seven poles 



ai = n{p + k,q2) 
"5 = ^{k,q2) , 



a2 = f2(p + k,qi + 92) , "3 = ^{p, 92) , 
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(D.25) 



Next we note that all double Born integrals except for can be found using combination of already 
calculated integral in Eq. (D.3) and the identity in Eq. (D.13). The results are following: 
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^The remaining three amplitudes in the figure 6 vanish. 
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The final integral we just work out: 
4 J 27r 27r 
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where in the first term of the remaining integral we have to close the contour in the upper complex 
plane, while in the second one we close above for 5z2 > Szi and below in the opposite case. Also, we know 
that 0^3 is in the positive complex plane, but we have to figure out the sign of imaginary part of (0^7 — 0^2) • 

Thus, for physical momenta p, fc, the second denominator in the last integral in Eq. (D.31) has a pole in 
the negative complex plane. With this in mind, we find the ^ for two cases of interest: 

lf\5z2 > Szl) = , (D.33) 
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Finally, taking the contact limit of the last integration in Eq. (D.31) we get: 

/f ^) = l[^\6z2 = 6zi + 0) = li^\6z2 = 6zi - 0) . (D.35) 

Next, in the contact limit and taking the average over medium, we have + q2± = 0. This makes some 
of the ai trivial. Directly from their definition we get: 

«2 = T , 0^4 = ^ = 0, 0^6 = 7^ = 0. (D.36) 

Po k+ 
Taking these equations into account we list all the contact limit integrals: 
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As discussed above, the seven inverse formation times reduced to only two combinations after averaging 
over the medium. We define these two inverse formation times in the following way: 

= «2 , = 0^7 + 0^3 . (D.43) 

Note that up until now we never used the soft gluon approximation, and all integrals and contact limits 
were exact. Finally, taking the soft gluon approximation gives: 

^ cjo 5 fis ^ cjo — oji . (D.44) 



E. T Wilson line 



(x"*", oo, oo) 



(x^, cx), oo) 



(x^, oo, oo) 



(x"*", oo, X_l) 





(x"*", oo, 



(x^, oo, x^) 



tI = P exp (^ig ^ dr U ^ 




Figure 12: The collinear and transverse gauge links. Left: the generic gauge, middle: regular covariant gauge and 
right: singular light-cone gauge. 



Wilson lines are a common ingredient for building gauge-invariant objects in field theory. For example, 
the collinear Wilson line Wn is used in SCET to dress the collinear quark field to form a gauge-invariant 
jet field Xn = Wn^n- However, in the light-cone gauge the collinear Wilson line disappears Wn = 1. In 
the context of SDIS it has been realized [60], [61] that additional transverse gauge link is required in this 
gauge in order to describe the final state interactions in a gauge invariant way. 

More recently, it has been proposed that SCET has to be expanded [54] by including a transverse 
gauge link for gauge invariance of certain non-perturbative matrix elements, like the transverse momentum- 
dependent parton distribution functions. Also, it has been shown in Ref. [54] that the SCET jet function 
calculated in the light-cone gauge in the presence of the transverse gauge link, which the authors denote as 
the transverse Wilson line Tn, is independent of the light-cone prescription, which unambiguously cancels 
between the gluon propagator and the gauge link, and leads to a result which is the same as in the covariant 
gauge. 

Xn^p-Tlwl^n^p, (E.l) 

Tl = Ve^p(^ig drA^-lj}j , (E.2) 

where V denotes the path ordering. As one can see from the figure 12 the effect of the transverse gauge 
link vanishes in the covariant gauge, since the gauge field is zero at infinity, while in the singular light-cone 
gauge this gauge link is non-zero depending on the boundary conditions, similar to the gluon propagator 
prescription in the light-cone gauge. 

Our calculation of the single and double Born amplitudes in the light-cone gauge in SCETq in the 
absence of the Wilson line leads to the same result as in the covariant gauge, but only in the —ie 
prescription. In all other prescriptions we get a different result, as one can see from Eq. (7.7) and table 7.2. 
In this appendix we show that once the Wilson line is introduced into the calculation, the prescription 
dependence cancels and the gauge invariance is restored similarly to the result in [54] for the SCET jet 
function calculated in the light-cone gauge. In the remainder of this section we assume the light-cone 
gauge, so Wn = 1 and we derive the Feynman rules for Tn Wilson line emission. 

A useful expression for the propagator of the gluon emitted from the transverse gauge link was derived 
in [54] and we show it in Figure 13, where the prescription dependence is encoded into the number C^^^^\ 



./i,6 



q 



^(Pres) 



q~^—ie 



Figure 13: The propagator of the T Wilson hne emitted gluon. 

Its dependence on the prescription is given in the table E. Note that the difference of this table and the 
propagator in figure 13 from [54] is due to the different notation of ingoing vs outgoing momentum flow 
into the vertex. 



Prescription 


\^] 


^(Pres) 




1 


1 


—is 


1 







q~^—i£ 




PV 


' \ ' ) 

2 \q++i£ ' q+-i£ J 


1 


ML 


1 

g++i£sign(g-) 





Table 3: Dependence of Coo^^^ on the light-cone prescription. 



First, from the form of the propagator in figure 13 one can see that the Wilson line cannot produce 
physical gluons in the final state, since h-e = in the light-cone gauge. However, this propagator contracted 
with the static source term i;^ doesn't vanish and is leading order in the effective theory power counting. 
In order to derive the Feynman rules of emission from the quark line we use the definition of the gauge 
invariant quark field Eq. (E.l) and the explicit expression for the transverse gauge link in Eq. (E.2). 
Finally we include the propagator in figure 13 to obtain first two Feynman rules in the figure 14 below. 

In order to derive similar Feynman rules with the gluon line we need the expression of the gauge 
invariant gluon field which is a straightforward generalization of the previous definition Eq. (4.19), now 
including the T Wilson line: 



^ 9 



TM(^dii + gA^)WnTj, (E.3) 



where the square brackets indicate that the derivative acts only within the brackets. From the expres- 
sion Eq. (E.3) and again using the definition Eq. (E.2) and the propagator from figure 13 we obtain the 
last two Feynman rules in figure 14. 
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Figure 14: Feynman rules for the gluon emission from the transverse Wilson line for single and double gluon 
emission from a quark and a gluon line. 
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